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With the continued improvement of sequencing technologies, the prospect of genome-based 
medicine is now at the forefront of scientific research. To reahze this potential, however, we 
need a revolutionary sequencing method for the cost-effective and rapid interrogation of individ- 
ual genomes. This capability is likely to be provided by a physical approach to probing DNA at 
the single nucleotide level. This is in sharp contrast to current techniques and instruments which 
probe, through chemical elongation, electrophoresis, and optical detection, length differences and 
terminating bases of strands of DNA. In this Colloquium we review several physical approaches to 
DNA detection that have the potential to deliver fast and low-cost sequencing. Center-fold to these 
approaches is the concept of nanochannels or nanopores which allow for the spatial confinement 
of DNA molecules. In addition to their possible impact in medicine and biology, the methods 
offer ideal test beds to study open scientific issues and challenges in the relatively unexplored 
area at the interface between solids, liquids, and biomolecules at the nanometer length scale. We 
emphasize the physics behind these methods and ideas, critically describe their advantages and 
drawbacks, and discuss future research opportunities in this field. 



Contents 



I. INTRODUCTION 



I. Introduction 
II. Current sequencing techniques 

III. Physical characteristics of DNA 

A. Structure 

B. Electronics 

IV. Nanopores and polynucleotides 

A. Nanopore-Polynucleotide concept 

B. Nanopore characteristics 

1. Biological pores 

2. Synthetic pores 

C. DNA translocation 

1 . Universal 

2. Specifics 

D. Pore and ion electronics 

1. Ionic and blockade currents 

2. Electrostatics 

V. Sequencing and detection 

A. Electronic 

1. Ionic blockade 

2. Transverse electronic current 

3. Capacitance 

B. Optics 

C. Force 

VI. Conclusions 

Acknowledgments 
References 



i 

i 
i 



i 

ia 
m 
m 

13 

M 

M 
M 
M 
M 

m 

[13 

m 
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sequencing improvements (jChanl . 120051 ) . genome-based 
medicine has come ever closer to reality. A low-cost 
and rapid method to sequence DNA would dramatically 
change the way we do medicine, and would give us new 
tools to study biological functions and evolution. 

This new technology would allow us to sequence and 
compare a plethora of individual human genomes, en- 
abling us to locate sequences which cause hereditary dis- 
eases and to discover mutated sequences. Individual 
medical consumers can then be tested for these known 
genetic defects. Thus, medicine and treatments can be 
tailored to their specific condition(s). This latter goal is 
known as personalized m,edicine and would be a tremen- 
dous advancement in the medical field. Furthermore, the 
genomic information gathered could be used in the study 
of biology by, e.g., comparing genomes between species 
to find common properties and functions. This informa- 
tion could be used both medically (e.g., in animal testing 
of potential cures for humans) and also to examine the 
evolutionary heritage of all species. 

In order to reach the goal of a rapid and low-cost se- 
quencing method, one cannot rely only on current tech- 
niques, which involve costs of about 10 million USD 
and sev eral months time to sequence a single human 
genome ( Fredla ke et all . l2006l ) . Improvements of current 
technology face both fundamental and practical limita- 
tions, such as a small - on the order of a 1000 base s 
- read- length limit using electrophoresis ( Chanl . |2005( ). 
which will restrict the impact of subsequent develop- 
ments in this direction. Therefore, radically novel ideas 
need to be implemented and demonstrated to be cost- 
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effective and accurate. 

There are several candidates which may 



fill this 



role (lAkeson et al. . 19991 : Deamer and Akeson 



2007 



Gracheva et al 



20001: iDeamer and Brantonl 1200? 
2006bl : iHeng et all 2005at Kasianowicz et al. 



Lindsav et 



19961 I2OOII: iLageravist aU 120061: iLee and Melle" 



20061: INHGRJ . 120061: 

Zwolak an d Di Vcntra, 12005 ). All these candidates 
have one trait in common: they employ nanoscale 
probes to examine the structural or electronic signatures 
of individual DNA bases. That is, they rely on physical 
differences between the bases. This is a major departure 
from existing sequencing paradigms which rely on 
chemical techniques and physical differences of strands 
of DNA. 

Most importantly, these proposals challenge our un- 
derstanding of, and ability to manipulate and probe, 
physical processes at the interface between solids, liq- 
uids, a nd biomolecules down to the nanometer scale 
regime (|Di Ventra et "all. I2004D . Indeed, in order to un- 
derstand the feasibility, speed, and accuracy of these 
novel approaches we are naturally led to examine sev- 
eral physical questions about the individual bases and 
the influence of the solid/liquid environment: 

• What is the difference in magnitude of physically 
measurable properties between the bases? 

• How do the nucleotide structural dynamics affect 
the measurable signals? 

• How do the different bases interact with the 
components of the detection apparatus, e.g., the 
nanopore, the surfaces, the electrodes, a scanning 
probe tip, or the other molecules present? 

• How does the atomic makeup and structure of the 
different bases affect the surrounding fluid and ionic 
dynamics? And vice versa, how do the latter affect 
the structure and electronics of the bases? 

• What are the significant sources of noise? 

In addition, many of the suggested sequencing meth- 
ods rely on nanopores either as a housing to contain the 
nanoscale probe(s) or as a restriction that causes differ- 
entiation in some signal between the bases. Thus, the 
fabrication of the nanopores and the DNA translocation 
dynamics have an important bearing on the following 
questions: 

• How well can one probe the DNA on the single base 
scale according to the dimensions and "uniformity" 
of nanopores? 

• How fast can one probe base differences? 

• What are the limits on read-lengths? 

• How fast and regular does the DNA translocate 
through the pore? 



The above points beget even more basic and general 
physical questions: 

• What is the meaning and role of electronic screen- 
ing at the nanometer scale and in a strongly con- 
fined and fluctuating environment? 

• What is the meaning of capacitance, thermal en- 
ergy, charging energy, etc., under these atypical 
conditions and at small length scales? How do these 
quantities evolve into their respective bulk proper- 
ties? 

• Do liquids show unexpected dynamical features at 

the nanoscale? 

• How do electrons move in "soft" materials and dy- 
namical environments? 

These questions will accompany us for the full length 
of this Colloquium. We stress their importance for the 
detection and sequencing approaches, review some par- 
tial answers found in the existing literature, and point 
out possible future research directions to explore them in 
more depth. 

The Colloquium is organized as follows: In Sec. |TT] 
we give a very brief account of current sequencing tech- 
niques. This primer will help the reader get familiar with 
the state of the art in this field. In Sec. IIIII we outline 
the physical properties of DNA and its bases. In Sec. lIVi 
we discuss nanopores as a useful building block for rapid 
DNA sequencing and detection. After this, we move on 
to physical approaches to DNA sequencing in Sec.|Vl Fi- 
nally, we conclude in Sec. |Vll 



II. CURRENT SEQUENCING TECHNIQUES 



A thorough introducti on to e xisting segucncing meth- 



ods can be found in 



io existing seqi 
iCha^. 120051 . rPredlake et all . 
I2OO6I . iDovichi and Zhangl l2000l and iJGll . 12004 refer- 
ences therein provide a more technical account. Present 
day sequencing methods are a n impr oved version of the 
Sanger method ( Sanger et all Il977( ). The sequencing 
proce ss can be divided into four overall steps (jChanl 
I2OO5I ): (1) DNA isolation, (2) sample preparation, (3) se- 
quence production, and (4) assembly and analysis. Step 
(1) is simply the isolation of the strand of DNA which 
needs to be sequenced. Within step (2) the DNA needs 
to be replicated and also broken into (many) very short 
strands. The length of the strands is dictated by the 
actual sequencing technology used. 

Step (3) combines three components for the detection 
of the bases in the DNA sequence, as shown in Figure [T] 
First, chemical elongation creates labeled strands of DNA 
with the random inser tion of a chain-ter minating nu- 
cleotide (introduced bv lSanger et all\l97'^ . Second, an 
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FIG. 1 (Color in online edition) The three processes used to 
produce the sequence of a strand of DNA. One starts with a 
DNA strand to be sequenced, called the template, to which 
one adds a primer which is complementary to part of the 
template. A DNA polymerase reaction extends the sequence 
starting from the primer by adding nucleotides to the 3' end. 
In the presence of a smaller amount of dideoxynucleotides, 
the chains terminate at various places along the complemen- 
tary strand, and after denaturation, produce single strands 
of different length. These strands are then sorted with cap- 
illary electrophoresis and detected by laser excitation of flu- 
orescent tags (one can either use multiple lanes each with 
its own dideoxynucleotide and fluorescently tagged primer, or 
use fluorescently tagged dideoxynucleotides). The detection 
information is then sent to a computer for assembly and post- 
processing. The terminology within this caption is defined in 
SecHLll 



fundamental barrier to the max imum read-l ength achiev- 
able when using electrophoresis ( Chanl . [2005h .^' The read- 
length is the longest strand of DNA which can be se- 
quenced accurately and efficiently within step (3). The 
read-length is limited because electrophoresis is sensi- 
tive to the physical difference between different lengths 
of (single-stranded) DNA. Thus, intuitively, one expects 
that as the strand gets longer, distinguishing a strand of 
N bases with one of -I- 1 bases becomes increasingly 
difficult because the percent difference in the strand prop- 
erties tends to zero. As a consequence, if a sequencing 
scheme allows for longer read-lengths of the DNA, this 
simplifies sample preparation and post-processing. Since 
the physical schemes described in this review rely on sin- 
gle nucleotide detection, the amount of sample prepara- 
tion should be considerably reduced. The electrophoresis 
step is also intrinsically slow ( Chanl . |2005| ) . 

This review focuses on methods which improve/modify 
step (3) of the sequencing process. However, as we just 
saw in the case of the read-length, the four sequencing 
steps are not independent. Modifications to step (3) can 
reduce the time and complexity of steps (2) and (4) and 
correspondingly reduce their costs. Thus, the techno- 
logical motivation for improving step (3) comes from its 
pivotal role in the sequencing process. 



III. PHYSICAL CHARACTERISTICS OF DNA 

In this section we give an introduction to the basic 
properties and structure of the polynucleotides (PN)* 
and their constituents. The bases are very similar to 
each other and, thus, these properties are crucial to un- 
derstanding their distinguishability via size, electronic 
states, or interactions with the surroundings. 



electrophoretic* process spatially separates the different 
lengths of DNA in a porous matrix. Third, an optical 
readout detects the fluorescent end groups or primers, 
which indicates the last base on each of the different 
lengths of DNA. 

Step (4) is the post-processing of the sequence data, 
which involves the reassembly of the short sequences to 
get the complete sequence of the original strand of DNA. 
Because of this, the short strands of DNA are required to 
have large overlapping sequences in order to match them 
up. 

From above, one sees that the current methods rely on 
very complex sample preparation and post-processing of 
the data. One of the main causes of this complexity is a 



A. Structure 

Both deoxyribonucleic acid (DNA) and ribonucleic 
acid (RNA) are built up of different bases attached to 
a sugar-phosphate backbone. The five bases that make 
up these polymers are shown in Fig. [51 These can be 
classified into two categories: the purine bases (A and 
G) and the pyrimidine bases (C, T, and U). The purine 
bases consist of a six- and a five-membered ring with a 
common edge. The pyrimidines have just a six-membered 
ring. These classifications are natural based on the chem- 
ical structure. We will see that for many physical prop- 
erties the distinction into purines and pyrimidines is not 
helpful. However, one would expect that some physi- 
cal properties would divide along this classification. For 



"Electrophoresis" is a general term meaning the action of driv- 
ing charged molecules/particles in a solution with an electric 
field. More specifically, it is taken to mean driving the parti- 
cles through a porous matrix. 



t Practical application of the chain-terminating chemistry can also 

limit read-lengths. 
* One also refers commonly to oligonucleotides, indicating a short 

strand of DNA/RNA. 
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FIG. 2 (Color in online edition) Atomistic structure of the 
nitrogenous bases found in DNA and RNA. For DNA, these 
are the purine bases adenine (A) and guanine (G), and the 
pyrimidine bases thymine (T) and cytosine (C). For most ri- 
bonucleic acids, uracil (U) takes the place of thymine. Each 
base is shown so that the bottom-most hydrogen indicates 
where the base attaches to the sugar group of the backbone. 
Also shown are the directions of their dipole moments. 



instance, purines are larger and thus one might expect 
that they would interact more strongly with surfaces in 
a confined space. 

The backbone structure of a polynucleotide is shown 
in Fig. [3^. Each monomer unit is a nucleotide, which 
consists of a base, phosphate group (PO4), and sugar. 
For DNA (RNA), the latter is the deoxyribose (ribose) 
sugar shown in Fig. [5}d(c). The difference of only a sin- 
gle hydroxyl group (OH) between RNA and DNA can be 
important for the global structure of the polynucleotide. 
Within the repeat unit of the polynucleotides there is 
one phosphate group, i.e., a nucleotide monophosphate. 
However, triphosphate monomers are used in chain ex- 
tension reactions. The formation of the double-stranded 
(ds) DNA helix occurs by the Watson-Crick base pairing, 
where A pairs with T, and G pairs with C. 

One of the most important properties of polynu- 
cleotides is that they are charged in solution. The pK^ 
of the phosphate group, i.e., the measure of how read- 
ily that group will give up a hydrogen cation (proton), 
is near one. Thus, under most ionic conditions (includ- 
ing physiological pH) , the backbone will contain a single 
negative charge for each nucleotide unit (or two nega- 
tive charges for a Watson-Crick pair of nucleotides in 
a double strand). In solution, though, there will be, 
on average, nearby counterions such as sodium (Na"*"), 
potassium (K"''), or magnesium (Mg^+), which neutral- 
ize a part of this charge. The fact that the strand of nu- 
cleotides is charged is what allows one to pull the DNA 
through nanopores with an electric field. 




FIG. 3 (Color in online edition) The atomistic structure of 
the nucleic acid backbones, (a) Two repeat units of the de- 
oxyribonucleic acid backbone. The phosphate groups have a 
negative charge in solution, which is shared between the two 
oxygens solely bound to the phosphorous atom. Also shown 
is the numbering system. The sugar carbon attached to the 
base is labeled 1' (the prime indicating that the numbering 
is on the sugar), then the other carbons are labeled around 
the ring up to 5'. This is the origin of the terminology "3'", 
"5'". A single strand can end after the 3' carbon or the 5' 
carbon. The rest of the figure gives the atomistic structure, 
as a unit in a polynucleotide, of the (b) deoxyribose sugar, 
(c) ribose sugar, and (d) dideoxyribose sugar. The latter is 
used in the Sanger method, see Sec. [Ill because the absence of 
the 3' hydroxyl group will not allow further chain extension. 
The typical nomenclature is to add a "d" in front of the base 
to indicate that it is attached to a deoxyribose sugar instead 
of a ribose sugar, e.g., poly(dA) is a polynucleotide of DNA, 
while poly(A) is a polynucleotide of RNA. 



There are other important properties that will help 
us understand the experiments and theoretical proposals 
below. For instance, a polynucleotide has a global ori- 
entation, with one end a 5' and the other end a 3', as 
described in Fig. [31 We will see in Sec. |TV] that this is 
important for the structural dynamics of polynucleotides 
translocating through a nanopore. 

There is also a global property called secondary struc- 
ture. For ds-DNA, for instance, there exist several differ- 
ent types of helices. The most common ones are called 
A-DNA and B-DNA. What changes the global structure 
of DNA between A- and B-DNA is the ionic and water en- 
vironment. B-DNA is preferred in aqueous environment 
because water molecules can bind in the groves along the 
helix. The base pair-base pair distance in the B-DNA 
helix is 3.4 A and there is a 36 degree angle between 
them, which gives about 10 bases per turn of the helix. 
The diameter of B-DNA is ~ 2 nm. A-DNA has different 
physical dimensions. The difference between B-DNA and 
A-DNA is probably not important for sequencing, but 
some proposals for detection of ds-DNA will be affected 
by such a change in global structure due to the environ- 
mental conditions. Furthermore, there is a process called 
denaturation where the two strands in ds-DNA unbind 
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FIG. 4 (Color in online edition) Schematic of a secondary 
structure of ss-RNA. On the top, stacked bases forming a 
single-stranded helix due to enthalpic effects (denoted as H 
in the figure). Entropic effects (denoted with S in the fig- 
ure) transform this structure into the one on the bottom, a 
random coil polymer. The homogeneity of the strand is im- 
portant for secondary structure, as are the ionic conditions 
and the temperature. For instance, poly(C) at room temper- 
ature and ne utral pH is mostly a single-stranded helix 1.3 nm 
in diameter (jArnott et aZ.l[r976l ). 



into single-stranded (ss) DNA molecules. 

Single strands show secondary structure as well, 
which also depends on ionic conditions and tempera- 
ture. A schematic of secondary structure in ss-RNA 
is shown in Fig. [?1 Generally, one can think of sec- 
ondary structure as a result of the competition be- 
tween enthalpic and entropic factors. The interac- 
tion energy of base stack ing is the main factor fa - 
voring secondary structure (jSearle and Williamsl . Il993t ) , 
such as the helix in Fig. [H On the other hand, the 
entropic factors are mai nly due to rotations of b ack- 
bone degrees of freedom (jSearle and Wimam£l . [T993h . If 
one raises the temperature, entropic factors will dom- 
inate and bring the strand to a random coil form. 
When this transition happens depends on the bases in 
the st rand (IVesnaver and Breslaueii Il99l[ ). ionic condi- 



tions (IDewev and Turne: 



19791) . and other environmen- 
tal factors (|Freier et al\ . 11981 ). The effect of secondary 
structure has alread y been seen i n the i onic currents 
through nanopores (|Akeson et ai\ . Il999l : iMeller et ai\ . 
[200Q) . 

Another type of secondary structure is a hairpin. 
RNA/DNA hairpins are single strands of RNA/DNA 
that wrap around to form a double strand. However, 
in doing so, they have to form a loop of unpaired bases, 
which causes an unfavorable strain on their formation. 
The double-stranded portion is called the stem, and 
sometimes a single-stranded portion can exist on the 
chain as well. This latter type of hairpin is used in some 
of the experiments described below. 

Let us also add that in some of the proposals for se- 
quencing, an exonuclease could be used to chop the DNA 
up into individual nucleotides. A nuclease is an enzyme 
that acts as a catalyst to polynucleotide breakup, e.g., 
the linking oxygen between nucleotide repeat units can 



TABLE I Sizes of the DNA bases (Vs) and nucleotides [Vn) 
in A^. Also given are the surface areas Ab and An in A^. 
We define the fraction of free volume left, Fd = {Vp — Vn) /Vp, 
where the pore volume, Vp — Trld^/A:, is for a cylinder of 
length I and diameter d. We assume each nucleotide occupies 
a length / ~ 7A, which is its own length in a completely ex- 
tended strand. We consider pore diameters of 15 and 20 A 
which are approximately equal to the a-hemolysin pore diam- 
eter in its stem. For comparison, the backbone and U volume 
(surface area) are 214 (225) and 128 (142), respectively. 

Base 



Vb {Ae 



Vn (A, 



Fi5 (Fa, 



A 


157 (166) 


349 (340) 


0.72 (0.841) 


G 


168 (177) 


359 (351) 


0.71 (0.837) 


C 


133 (147) 


324 (319) 


0.74 (0.853) 


T 


150 (163) 


339 (331) 


0.73 (0.846) 



be hydrolyzed to two OH groups. An exonuclease is a nu- 
clease that "eats away" at the ends of the polynucleotide, 
breaking off one nucleotide at a time. 

Some particular properties of interest to us are the 
base and nucleotide sizes. The sizes of the DNA bases 
and corresponding deoxyribonuc leic acids are shown in 
Table HI Using geometries from IZwolak and Di Ventral . 
120051 . these sizes were computed by drawing a sphere 
around each atom of its van der Waals radius, and tak- 
ing the union of these spheres. From the table, one can 
observe that even in the most ideal case, the size of the 
base alone is unlikely to provide a distinguishable signal. 
This will be discussed at length in Sec. IIV.DI 



B. Electronics 

Two detection/sequencing schemes propose to use 
transport and voltage fluctuations to distinguish the 
bases. These will be sensitive to the electronic struc- 
ture of the different nucleotides and also to the geometry 
and local environment. 

Transport through nucleotides will detect differences 
in electronic states via their energy and spatial exten- 
sion. To understand the differences between the bases we 
compute the density of states of the nucleotides (placed 
between two electrodes) projected onto the backbone and 
bases, as shown in Fig. \5[ The plot shows that the dif- 
ferent nucleotides do not have a considerably different 
electronic structure, having molecular states very near 
in energy compared to the position of the Fermi level. 
Most of the density of states around the Fermi level is 
contributed by the backbone. This is in part due to the 
upright geometry in Fig. Eljb), which allows for a large 
contact area of the backbone with an electrode. The 
different spatial extension of the bases causes their con- 
tribution to the DOS at the Fermi level to differ. These 
differences will be heavily influenced by the nucleotide 
geometry and orientation. 

However, the basic input to determine the base- 
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LUMO 



jbi- * * * * 




FIG. 5 (Color in online edition) (a) Projected density of 
states of the DNA nucleotides between two gold elec t rodes 
shown in (b), as computed in IZwolak and Di Ventral . |2005| 
with tight-binding parameters. Ef is the Fermi level of gold 
within the same approach. 



TABLE II Theoretically calculated dipole moments of the 
bases (pb) and the deoxyribonucleic acids (pjv) in Debye 
0.21eA). For comparison, the dipole moment of U is 
4.79 D (p^'' = 4.2 D) and the backbone dipole by itself is 
1.97 D. Experimental values, p ^'', are from Kulako ws ef aU . 
11974 IWeber and Cravenl . ll990l . and lDevoe and Tinoc"cl . ll962l . 



Base 



Pb 



exp 

Pb 



Pn 



A 
G 
C 
T 



2.33 
7.17 
7.22 
4.72 



2.5 
7.1 
7.0 
4.1 



4.76 
7.76 
8.55 
7.56 



electrode coupling is the character of the molecular 
states. For the isolated bases, this is shown in Fig. [H] 
for the highest-occupied molecular orbital (HOMO) and 
lowest-unoccupied molecular orbital (LUMO). For nu- 
cleotides only weakly contacted with electrodes, the char- 
acter of these states is going to determine how well the 
nucleotide can couple to charge carriers in the electrodes. 
We can see that for all the bases, the states are dis- 
tributed around the rings. These states remain that way 
even in the presence of the passivated or charged back- 
bone (although they may not remain the HOMO and 
LUMO states of the nucleotide). The spatial extension 
of the base is roughly correlated with the relevant elec- 
tronic states and, thus, will determine a large part of the 
couphng of the nucleotide to the electrodes. 

On the other hand, other electronic approaches may 
measure the dipole and higher moments. The dipole mo- 
ment magnitudes are shown in Table |TT1 with the corre- 
sponding directions shown in Fig. [21 The dipole moments 
are for the isolated DNA bases and for the correspond- 
ing nucleotides. The moments and the molecular wave- 
functions wer e computed within Hartree-F ock using the 
geometries in IZwolak and Di Ventral [200I 

Finally, a very important property of PNs relevant to 




T 



FIG. 6 (Color in online edition) Isosurfaces of the HOMO and 
LUMO states of the isolated bases. The red and blue colors 
indicate opposite signs of the wavefunction. The HOMO and 
LUMO states do not change from the individual base to the 
passivated nucleotide. However, when the backbone becomes 
charged in solution, the HOMO and LUMO states will shift 
in both character and energy. 



the studies discussed in this review, relates to its bond- 
ing properties at surfaces, and in particular at the inte- 
rior surfaces of nanopores. These interactions, and their 
effect on the PN dynamics in confined geometries, is a 
property of both the nucleotides and the type of surface, 
thus making it a very complex issue. We defer its discus- 
sion to later sections on particular experiments. 



IV. NANOPORES AND POLYNUCLEOTIDES 

About a decade ago, Kasianowicz and collabora- 
tors were able to pull single-stranded polynucleotides 
through a biological nanopore by applying a voltage 
across the por e which p u lls on the charged PN back- 
bone (Kas ianowicz et al\ . Il996( ). These authors have 
detected the translocation of the PN via measurement 
of the blockade current, to be discussed below. Since 
this pioneering experiment, nanopores have been used 
to extract a variety of information characterizing the 
translocation, dynamics, and interactions of both single- 
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and double-stranded PNs in nanopores (lAkeson et al. 


199g 


ilButler et al.. 2006; Chang et all 2004; Chen et al. 


2004 


: Deamer and Brantonl 12002; Foloeea et al. 1 


005a 


Henrickson et all 2000|: 


Kasianowicz et all 


1996 


2001; Mathe ei aL. 20051 


Meller et al. 2000l 


2001 


Storm et all. 2005bl). Further progress is heme made 



by moving to synthetic nanopores where techniques are 
being developed to c ontrol nanopor e shapes, sizes, and 



other characteristics (iLi et al. UgOOal I2OO1I : IStorm et all . 
I2003t IWanunu and Melleil . l2007^ . Successful sequencing 
and detection will require control of these characteristics. 

In addition to experiments, theoretical and com- 
putational work is underway to help understand 
polynucleotide translocation through nanoscale pores. 
There are two fruitful approaches to this complex issue: 
phenomenological models and molecular dynamics. The 
phenomenological models provide a highly reduced 
description of the polymer dynamics, but they are able 
to elucidate the dependence of the dynamics on param- 
eters such as t he polymer length, p ore di r nensions, and 
appli ed field ([Ambiornsson et aZ. I2OO2I; 'Chern et al), 

Kong and M uthuku mar , 
'l99S; 



2001 



2002 



Chuang et al 



Luo et al. 



Loebl et al 
200^ 



20021: 

20031 I Lubenskv and Nelsonl . 



Matvsiak et all 120061: iMelleii, 



2003 ; 



Muthukumaij. Il999l 120011: iSlonkina and Kolomeiskv . 
2003t ISung and Parkl Tl996D . On the other hand, if 



one wants to understand how to probe physical differ- 
ences between the bases, an atomistic description of 
the polynucleotide dynamics is necessary. Molecular 
dynamics simulations coupled to other computational 
methods have been used in this context to study the 
signals and fluctuations expected when mcasviring 
different physical quanti ties as the PN translocates 
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In this section, we first review the basic concepts of the 
nanopore-polynucleotide experiment. We will then move 
on to discuss some experimental results on biological and 
synthetic pores. We also include a number of interesting 
physical results on PN translocation and pore electron- 
ics, which demonstrate the wealth of fundamental physics 
that is contained in this new field. 



A. Nanopore-Polynucleotide concept 



Prior to [Kasianowicz et al. 1. 119961 . researchers were al- 
ready interested in the ability to pull small charged 
mole c ules and polymers through ion channels (iBavlev , 
1994'; Bezrukov and Kasianowicz, '1993'; Bezruko v et al . 
"996, 1994; Bustamante et al, 1995; Henry et a/.i , ll98£ ; 
Kasianowicz and Bezrukovl 1 1 995l ) . Thus, the basic ques- 
tion was whether PNs could also translocate through, 
and be detected by, a pore. A schematic of this process 
is shown in Figure [71 The figure shows a membrane/thin 
layer which divides a solution into two halves. An elec- 




FIG. 7 (Color in online edition) Schematic of the nanopore- 
PN experiment. A membrane divides the solution into what 
is commonly called the cis and trans chambers. A bias across 
the separation creates a field across the membrane and 
drives an ionic current through the pore as shown in (a). 
The field also pulls on the negatively charged PN backbone, 
which causes the PN to get captured by and then translocate 
through the pore as shown in (b) . While the PN is within the 
pore, ions are partially prevented from occupying and flowing 
through the pore, thus reducing the ionic current. The pore 
is characterized by its length Lp and average diameter Dp. 



trode is placed into each halve. Without the PN, the 
electrodes pull the ions through the channel to create an 
ionic current Iq of the open channel (see Figure [8|). 

Due to the negative charges on the phosphate groups, 
the PN is pulled to the positively biased (cathode) half 
of the solution. Eventually the strand is captured and 
enters the pore. During a time td, the translocation du- 
ration, the strand partially blocks ions from the pore as 
shown in Figure W[i>). Even though the nucleotide is 
charged, it carries very little ionic current through the 
pore because of its slow velocity compared to the ions. 
Thus, its presence gives the blockade current lb shown 
in Figure [H so that the PN translocation event can be 
detected. The translocation events also provide a way to 
obtain information about the PN, such as its length and 
limited information on composition and dynamics. 

The actual values of the ionic current of the open pore 
and the blockade current are due to several factors. The 
open pore resistance is determined by both an energy 
penalty (mainly electrostatic) and an entropic barrier to 
bring charges into the pore. We will see below in sub- 
section IIV.DI that this is quite a complicated problem in 
itself, due to the microscopic details of electric fields in 
reduced geometries and fiuctuating environments (e.g., 
the value of a screened charge in an electric field at the 
nanoscopic level can not be simply viewed as a charge 
within a dielectric environment), and also the surface 
physics of pores in solution. The blockade current is thus 
also a complex phenomenon. We will discuss this issue 
in the context of the particular experiments described 
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FIG. 8 Example of an ionic current signal. Upper Figure: 
A voltage of 120 mV drives a current through a a-hemolysin 
pore. When poly(U) is added to the cis chamber, it gener- 
ates blockade events. Lower Figure: Two examples of block- 
ade events. The average open pore current, Jo, the average 
blockade current, /;,, and the translocati on duration, tg, are 
all in dicated on the figure. Adapted from iKasianowicz et all 
Il996l . 
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FIG. 9 Slab view of the a— hemolysin pore. The uniform 
light gray area represents the lipid bilayer into which the 
Of— hemolysin assembles to make the pore. The initial part 
of the pore with a much wider diameter is the vestibule. 
The approximately 5 nm long neck is the pore stem. Some 
other characteristics of the pore are a ~18 nm^ stem vol- 
ume (jPeamer and Brantonl . [20021 ) . and at 120 mV, the pore 
current is 120 pA, giving almost 10® ions/s. 



below. 



1. Biological pores 



B. Nanopore characteristics 

The current state of the art is divided into two Hnes 
of research. One l ine uses biqloKic al pores, generally 
a— hemolysin pores ( Song et aLl . ll996t ). These pores have 
the right size scale to detect differences in strands of PN. 
However, tailoring characteristics, such as the pore diam- 
eter, is not straightforward. The other line of research 
operates with synthetic pores, which can be reasonably 
controlled down to the sub-nanometer range. One can 
also imagine making hybrid devices to help control pore 
properties. 

The characteristics of the nanopore and the ability to 
control them will be extremely important for sequenc- 
ing methods. For instance, methods which propose the 
use of nanoscale probes embedded in the pore will re- 
quire a pore size that maximizes the signal difference be- 
tween the bases beyond the many unavoidable sources of 
noise. In most cases this means a pore diameter with 
the same width as a single strand of DNA. However, the 
maximization of the signal difference has to be balanced 
with other effects such as DNA-surface interaction, which 
is minimal with a large diameter pore, and DNA cap- 
ture/translocation, which will not occur at small pore 
sizes. In addition, the pore has to be fabricated in a way 
to make possible the embedding of a nanoscale probe. 
This will put restrictions on the types of material for the 
pore and its sizes/shapes. 



The a— hemolysin pore is shown in FigurelHl with some 
characteristics given in the figure caption. In this bio- 
logical pore the smallest restriction that the PN has to 
translocate through is ^--^1.4 nm. Thus, ds-DNA, at a di- 
ameter of 2 nm, can not translocate throu gh the pore, 
but ss-DNA can ( Kasianowicz et a/] . Il996l ). The pore 
is also both small and long enough that the PN has to 
be locally extended, which gives an entropic barrier for 
transport due to unraveling of the polymer. We give an 
overview of some a-hemolysin experiments in Table IIIII 

A typical experimental setup with an a— hemolysin 
pore is shown in Figure 1101 There are two chambers, 
the cis and the trans, with a buffered solution, e.g., 
of KCl. Between the two is a Teflon partition with a 
small orifice where the lipid bilayer is formed. When the 
hemolysin subunits are added they spontaneously form 
into the pore. The formation of the pore is detected 
by the appearance of an ionic current between the two 
chambers. 

The experiment is then conducted by adding PNs into 
the cis c hamber. The first experim ents were done with 
poly(U) (jKasianowicz et aJ\ . Il996t) . Upon its addition, 
transient blockades of ionic current were observed. The 
blockade events actually fall under three different types 
distinguished by their lifetime, as shown in Figure [TT] 
There are fast blockades that are independent of the 
poly(U) length. Thus, most likely, these are events where 
the strands just cover the entrance to the pore or just 
partially enter the pore, but do not translocate through 
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TABLE III Overview of some a- hemolysin experiments. T represents the temperature and vdna is the DNA velocity through 
the pore. "p(X)" stands for a polymer of X, and L is the length of the polymer in nucleotides (nt). The designation 3' ^ 5' 
indicates an experiment which looked at the translocation direction. 

Ref. PN L (nt) Result 



Kasianowicz et al. . 1996 
Akeson et al. 1999 
Meller et al. 2000 
Meller et al . 2001 
Mathe et al. 2005 
Butler et al. 2006 



p(U), others ~ 150 - 450 
p(A),p(C),p(dC),p(U),p(A)(C) ~ 100 - 200 
p(dA),p(dC),p(dA)(dC),p(dAdC), others 100 
p(dA) 4-100 
p( A), hairpin ~ 50 

p(A),p(C),p(A)(C),p(C)(A) 50,75 



tdoc L 
discrimination 
td <x discrimination 
Ji, transition at L ~ Lp,VDNA 
3' ^5' 
3' 5', discrimination 




FIG. 10 A typical experimental setup for an a-hemolysin 
pore. The positive voltage is applied to the trans chamber 
and the negative voltage to the cis one. From lAk eson et all 
Il999l . 



it. However, surprisingly, both the other two types of 
events have lifetimes linearly dependent on the length 
of po ly(U), and inversely depe ndent on the applied volt- 
age. iKasianowicz eldl ■ Il996l conjectured that this could 
be due to different translocation speeds of strands enter- 
ing with either the 3' or 5' end. Later experiments con- 
firmed this conjecture by examining strands composed of 
two homogeneous blocks and exploiting the pore's abil- 
ity to distinguish between different nucleotides, see be- 
low (Mathe et al . .2005.) . 



2. Synthetic pores 

Synthetic pores offer additional opportunities for de- 
tection and sequencing of PNs. For instance, one can 
adjust the pore dimensions and properties to meet the 
needs of a particular experiment. Also, they open the 
possibility of integration of external sensors and probes, 
such as transverse electrodes. In addition, the param- 
eter range for their operation is larger (although the a- 
hemolysin pore is quite robust, biological pores in general 
are only open under certain voltages, ionic concentra- 
tions, and temperatures (Hillc, 2001)), and indeed differ- 
ent parameters such as salt concentration, temperature, 
voltage, and viscosity may be key to operating the pore 
in a regime where sequencing and detection is possible or 
most optimal, because these conditions control proper- 
ties like the translocation velocity and capture rate (see. 
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FIG. 11 Translocation dur ation (lifetime) of t h e ioni c block- 
ade events measured by IKasianowicz et "al. I 119961 . The 
event s fall under three diffe rent lifetimes, see the text. 
From IKasianowicz et al 1. 119961 . 



e.g.. iFologea et all [2005bl : iHenrickson et all l2000t) . We 

give an overview of some of the synthetic pores made to 
date and the translocation experiments in Table [TVl 

The fabrication of synthetic nanopores is still in its 
nascent stages. Mainly two groups have pioneered tech- 
niques for solid-state pore fabrication. Golovchenko, 
working with J. Li, has developed a technique using 
low-energy ion beam to s culpt ("ion-be am sculpting") 
a nanoscale hole in Si3N4 ( Li et a/.| , l200lh . Dekker et al. 
have developed a technique base d on an elec t ron b eam 
and a visual feedback procedure ( Storm et all . l2003l ) . 

The technique de veloped by Go lovchenko et al is 
shown in Figure fT2l ( Li et all 1200 ll ). One starts by cre- 
ating a large diameter pore in a solid-state membrane 
using a focused ion beam (Figure IT^). In their case, 
a ~ 60 nm pore in a Si3N4 membrane is first created. 
Then this pore is exposed to an Ar+ beam, which, in- 
stead of knocking atoms off the membrane and opening 
the pore further, ac tivates a diffus ion process and the 
pore starts to close ( Li et a/.l . l200ll ). The current of the 
Ar"*" coming through the pore is directly dependent on 
the pore area. Thus, one can measure the Ar"*" current, 
shown in Figure 112b .. as a function of time, and use this 
to controUably shrink the pore down to the nanometer 
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TABLE IV Overview of some synthetic nanopores. Next to each pore type, their diameters and lengths are given as {{Dp};{Lp}) 
in units of nanometers. 



Ref. 



Pore 



PN 



L (nt) 



Result 



Li et al. 2003 
Chen et al. 2004 
Storm et al. 2005b 
Fologea et al . 2005a 
Chang et al. 2004 



Si3N4 (3,10;5-10) 
SisNi (15;n/a) 
Si02 (10;20) 
Si3N4 (4;5-10) 

Si02 (4.4;50-60) 



Heng et al. 2004 Si3N4 (1,2.4;10,30) 



bio-ds-DNA 
A-DNA, other 
bio-ds-DNA 
bio- ss- and ds-DNA 
bio-ds-DNA 
p(dT), ds-DNA 



3,10 kbp 
8.5, 10, 3 kbp 
~6-97 kbp 
3 kbp 
200 bp 
50-1500 bp 



Folded DNA h events 
vdna oc V, folded 
tct oc 
denaturation 
h > lo 
length discrimination 




FIG. 12 Fabrication technique of iLi et oil 1200]] for creating 
nanoscale solid-state pores: a Si3N4 membrane with a large 
pore is first created (panel b). An ion beam is then focused 
at this large pore, activating a diffusion process which closes 
the hole (panel c). The current of Ar^ decreases as the pore 
shrinks (panel a), the monitoring of which can be used to 
control pore size. From lLi et al I I2OOII . 



scale (FigurefT^V The precise composit ion of silicon an d 
nitrogen around the pore is not known ( Li et aLl . l200ll ). 

Dekker and co-workers have developed a similar tech- 
nique based on the idea of pore diameter reduction upon 
irradiation ( Storm et al I I2OO3D . They use a high-energy 
electron beam and Si02 membrane. The initial, larger 
pore can be formed in different ways. But once the initial 
pore is made, a transmission electron microscope (TEM) 
is used to shrink it. By using the imaging capability of 



the T EM, one can visua lly watch the pore diameter re- 



duce (jStorm et ali . \200§t) . Since the rate of size reduction 
is very slow (0.3 nm per minute in these experiments), 
the visual images can be used to monitor the diameter 
of the pore and the process stopped when the desired 
size is reached. The resolution of the microscope (~0.2 
nm) then sets the limit of accuracy for reaching the de- 
sired size. However, due to the roughness of the surface. 



one c an only control the pore size to about 1 nm dimen- 
sions (jStorm et aU 120031 ). 

Synthetic pores have properties that differ from the 
biological pores. The surfaces can be charged in solu- 
tion, creating an additional complication to understand- 
ing ionic and PN transport. For instance, silicon oxide 
surfaces can have a negative surface charge density on 
the o rder of 10~^ e/A in aqueous solution ( Chang et all . 
I2004D . This would give quite a substantial charge within 
any reasonably sized nanopore, and would have to be 
neutralized by counterions. Biological pores, though, can 
have unusual potential pro files w hich may have internal 
sites with trapped charges ( Hilld . l2001i) . 

The observation that high energy beams can cause 
pores to shrink is a very interesting phenomenon in and 
of itself, in addition to its implications for the fabrica- 
tion of n anoscale structures/ pores. We want to mention 
here that Stor m et "all. 120031 have observed that there is 
a transition from shrinking to growing as the pore diam- 
eter is increased, and have explained this observation by 
a surface tension effect. 



C. DNA translocation 

One can divide the translocation of PN through a pore 
into two categories: universal properties of polymer dy- 
namics (entropic forces, Brownian motion, charges and 
screening), and specific properties that rely on the atomic 
compositions of the nucleotides (e.g., interaction poten- 
tials with the pore surface). Depending on the quantity 
under consideration, and the parameter regimes of the 
device, either one or both of these categories will be im- 
portant. 



1. Universal 

The two basic properties of charged polymer dynamics 
are the processes of capture and translocation. The cap- 
ture of the polymer will depend on the diffusion of the 
polymer from the bulk to the pore and on local effects 
around the pore, such as the electric field and interac- 
tions between the entrance of the pore and the polymer. 
The captu re rate will depend on concentration and a' 



plied bias ( Henrickson et al\ . I2OOOI : iNakane et al 
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as well as what molecule is under investigation, and will 
have repercussions on the ability to detect and sequence. 
The translocation through the pore will be driven by the 
applied bias, but depends on many factors, including the 
polymer-pore interactions, ionic effects, and viscous drag. 
However, two properties are common to translocation in 
nanopores, namely the effective charge and screening of 
polymers within the pore, and the related issue of where 
the applied voltage drops. 

In the absence of a polymer within a pore, one expects 
that the majority of the voltage drop between the two 
ion chambers occurs in the pore since it has a resistance 
higher than the surrounding solution. However, depend- 
ing on its shape (which depends on whether the latter is 
biological or synthetic), one would expect that the pres- 
ence of a polymer in the pore can change significantly 
the voltage drop of the system. For instance, with the a- 
hemolysin pore in Figure [9l one might expect that when 
the polymer creates a much higher resistance, this oc- 
curs primarily in the region of the pore which is roughly 
2 nm in diameter. Thus, the potential profile will change 
and drop more significantly over the stem. This would 
cause the applied bias to act most significantly on the 
nucleotides (or charged polymer units) that are located 
within the stem. 

There is some in direct evidence that this is the case. 
iMeller et all l200ll showed that there is a transition in 
velocity dependence of the translocation process. Above 
about 12 nucleotides, the velocity of the strands is in- 
dependent of their length, suggesting that the bias is 
pulling only on a finite number of charges and is acting 
solely against the drag of the polymer within the stem. 
Below about 12 nucleotides, the velocity increases (in a 
highly non-linear way) with decreasing length. For 12 nu- 
cleotides, the contour length is approximately 5 nm, thus 
in agreement with the pore stem length. This transition 
is explained by a model of polymer translocation throug h 
a pore of finite length i Slonkina and Kolomeiskvl |2003( ). 
which gives a transition in velocity when the polymer's 
length is equal to the pore length. 

Another interesting fact is that there is a transition 
in polymer dynamics based on the pore size and pore- 
polymer interaction. For pores with diameter roughly 
the size of the polymer width, one expects strong inter- 
actions of the polymer with the walls of the pore, thus 
strongly increasing its resistance to fiow. When this is the 
case, the timescale of polymer translocation, td, would be 
controlled by this resistance, which would dominate over 
other factors, such as polymer unraveling and drag out- 
side the pore (see Figure [T5|. In this regime, one expects 
the translocation duration to be linearly proportional to 
polymer length (for polymer lengths much larger than 
the po re length). Indeed this is what has been mea- 
sured ( Kasianowicz et aLl . ll996l ). On the other hand, if 
the pore is very wide (or there is very little polymer-pore 
interaction), the translocation time can be controlled by 
other factors, and one does not expect this time to have 
a linear relation with the polymer length. 




P 



FIG. 13 Drag of a co iled-up polymer outside a pore. 
From IStorm et aH l2005al . 



Such a nonlinearity h as been observed and ex plained 
by IStorm et all l2005al and IStorm et all l2005bi These 
authors observe that the (most probable) translocation 
duration scales as 



tdcc L 



1.27 



(1) 



for unfolded ds-DNA translocation. By considering a 
model that accounts for the driving electric force and 
the hydrodynam ic drag force on the coiled-up polymer 
outside the pore, IStorm et al. I. I2005bl obtain a transloca- 
tion duration that scales with polymer length as L^-^^, 
which compares well with the experimental value. The 
unusual exponents come from the fact that the radius 
(of gyration) of a polymer coil scales as function of its 
length with non-integer exponent, and depends on the 
polymer type (self-avoiding, etc.). The drag force will 
be dependent on the surface area of the coiled polymer. 
Most sequencing technologies will probably have to deal 
with small pore diameters, and thus in the regime where 
drag of the polynucleotides within the pore is significant. 
Therefore, in these devices the translocation duration 
should scale with the length of the polynucleotide. More 
generally though, this result shows that if one wants to 
use a nanopore as a polymer detector, the latter has to 
be first calibrated to take into account specific details 
which affect the measured characteristic signals (similar 
cali bration has to be perfo rmed in the case of sequenc- 
ing ( Lagergvist et aLl . l2006[ ). see subsection IV. A.2p . 

The linearly extended polymer may not be the only 
one able to translocate. If the pore diameter Dp is large 
enough to accommodate multiple strands, there could be 
ionic blockade events corresponding to folded polymers. 
In this case, let us analyze the forces on the DNA nu- 
cleotide pairs, and, for folded DNA, on the nucleotides 
that fall within the pore. This analysis highlights the 
complexity of polymer dynamics in nanopores. If the ap- 
plied voltage drops solely over the pore, the driving force 
on the polymer in a given region of the pore will be 



-^driving (^) ~ ^^-^dr 



(2) 



where m is the number of folds in the region, and ^X-iumg 
is the driving force on a single nucleotide pair. 



driving ^^ff 



E 



(3) 
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where Zeff « — 0.5e (see Sec. lV.C() is the effective charge 
on the nucleotide pair and E is the electric field. If pore- 
polymer interactions are responsible for the majority of 
the drag force, then one also obtains a similar relation 



Fdrag{m) « mf^-^g 



(4) 



where, e.g., F^^.^^ « r]v, with 77 a coefhcient that is pro- 
portional to the surface area for contact of the polymer 
and pore surface (that is the reason why 77 — > mrj for the 
folded polymer), and v is the polymer velocity. For a 
constant polymer velocity*. 



Fdrag — ^driving 



(5) 



and, thus, the factors of m will cancel, and v would be 
the same regardless of whether the polymer is folded or 
unfolded (v = Vo)- Now let us consider a polymer that 
has a region of length Lm of m folds, with ^ Lp. If 
we suppose volume exclusion to be responsible for ionic 
blockade (see next section), the blockade current of the 
folded and unfolded polymer will be related by 



lo - ir = - lb ) 



(6) 



where /q is the open pore current, is the blockade cur- 
rent for an unfolded polymer, and /™ is the correspond- 
ing quantity for the polymer with m folds. Given Eq. [51 
one finds that the total charge blocked for the translo- 
cation of the folded polymer region is equal to the total 
charge blocked by a linearly extended region of length 
mLm, which contains the same number of nucleotides, 



(/o-O — = (/o-4°) — 

V Vo 



(J) 



This is what ha s been dubbed "the ru le of constant event 
charge deficit" ( Fologea et al\ . [2005al ) 



ecd = A/fctrf = constant. 



(8) 



where A/b = lo — It, and td is the translocation time. In 
experiments where this holds, it would seem to justify an 
excluded volume model for the blockade current. 

This phenomenon has indeed been observed in exper- 
iments with 10 nm di ameter Si 3 N4 p ore (5-10 nm thick) 
and 10 kbp ds-DNA (|Li al\ . l2003t) . A smaller Si3N4 
pore of 4 nm diameter also showed a set of events that 
were d istributed according to this "rule" I Fologea et all . 
l2005ah . However, these observations were only for larger 
pores (4 and 10 nm for ds-DNA, 4 nm for ss-DNA), and 
thus they are not necessarily indicative of the physics 
of smaller pores. There ar e also events that do not fall 
on curves of constant ecd ijFologea et IIl,|2Q05i)- These 



This, of course, assumes that there is such a thing as a "constant 
velocity" for this nanoscale object. This is not strictly true, and 
estimates of drag and subsequent analysis may be affected by a 
varying velocity. 



seem to be events where the ds-DNA temporarily bond to 
the walls of the pore or have some other factor controlling 
their velocity (besides viscous forces) , such as untangling 
of the ds-DNA outside of the pore. 

To get the "constant event charge deficit," one has to 
assume that a pore-polymer interaction or, at least, a 
polymer interaction specific to the pore region (e.g., an 
electrostatic drag induced by the confinement) is caus- 
ing the drag which will give the relation in Eq. [4l With 
this type of drag, one expects the translocation dura- 
tion to scale as L, the length of the polymer. How- 
ever, IStorm et all l2005al find that the (most probable) 
translocation duration scales as in Eq. [T] This seems to 
indicate that it is not the polymer-pore interaction creat- 
ing the drag but rather the externally jumbled polymer. 
There are many possible physical explanations for such 
a discrepancy and other experiments are needed to have 
a complete understanding of these effects. 



2. Specifics 

Here we give just one example of an interesting 
specific property observed for PN experiments with 
Q- hemolysin pores. We previous ly discussed results 
of lKasianowicz et al. I, Il996l (see also lAkeson et al\ . Il999t ) 
which raise an intriguing question about the direc- 
tionality of translocation and capture of the PN. Are 
these two processes dependent on which end, the 5' or 
3', translocates first? Recently, both experiments and 
MD simulations have confirmed direc tionality-dependent 
translocation. iMathe et all l2005l performed experi- 
men ts on DNA hairpins, a s shown in Figure [T4| (see 
also iHenrickson et all i2000l ) . The experiment has been 
performed by pulling the hairpin into the pore, measuring 
the ionic blockade, then turning off the pulling voltage for 
a time toff and using a small probing voltage to detect 
whether the pore was still blocked. These authors have 
found that the hairpins with the 3' end dangling (translo- 
cating 3' 5') has a lower blockade current than the 5' 
end (translocating 5' 3'), see Figure [Ml Further, these 
two blockade currents are in good agreement with the two 
peaks found from performing the typical single-stranded 
experiment (with no hairpin). This indicates that early 
experiments obtaining translocation events with two dif- 
ferent blockade currents were measuring the difference in 
orientation of the strands. In addition, the diffusion and 
velocity (when pulled) of the hairpin out of the pore was 
found to be much slower for the 3' end entering, e.g., 
when the strand is translocating 5' — > 3'. This also is a 
good indicator that the early observation of two lifetimes 
(with different capture rates) of translocation events was 
likely due to different orientation of the translocating 
DNA. 

In addition, iMathe et all i2005i performed MD sim- 
ulations which demonstrated that in a pore the bases 
on the DNA tilt toward the 5' end. This explains the 
slower motion of the DNA, because when moving 5' end 
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FIG. 14 Schematic of a hairpin experiment with an a- 
hemolysin pore, (a) A typical experiment on ss-DNA can 
pull the strand through from 3' — > 5' or from 5' 3' leading 
to blockades with different characteristics, (b) The double- 
stranded portion of the DNA hairpin can not transverse the 
pore. Thus, the hairpin can be pulled as far as possible 
into the pore, and one can turn off the pulling voltage for 
some time t^f / after which one can check whether the strand 
is stil l in the pore by ap plying some small probing voltage. 



IS stn l m tne pore by ap 
From iMathe et all 120051 . 



first there is additional "mechanical" friction. This phe- 
nomenon will be important not only for detection of dif- 
ferent types of polynucleotides, but also for any of the 
sequencing proposals described below. For instance, if 
there is a different average base orientation depending 
on which end enters first, the current distributions of 
the bases (see subsection IV.A.2|) may also be different 
depending on whether transverse control erases the di- 
rectional differences. One may also ask whether this dif- 
ference in 5' and 3' orientations will be present in a syn- 
thetic pore, or whether it depends on properties specific 
to the a-hemolysin pore. Further, the interplay between 
this effect and that of secondary structure is still unclear, 
pointing to the need for more experimental and theoret- 
ical investigation. 



D. Pore and ion electronics 

From the above discussion we conclude that ionic cur- 
rents through the nanopore are important for the detec- 
tion of polynucleotides. However, much is still not un- 
derstood about the translocation of ions themselves. The 
interplay between volume exclusion, steric blockage, elec- 
trostatics, solution-ion dynamics, and interaction with 
the pore itself is quite complex and poorly understood. 



1. Ionic and blockade currents 

In part, the blockade current forms due to "volume ex- 
clusion" . The polynucleotide occupies part of the pore's 
volume which partially blocks ions from occupying the 
pore. This reduces the number of charge carriers and 
thus the current. How much of the volume is occupied 



by the PN depends on the structure and composition of 
the strand. PNs with helical secondary structure could 
block more of the pores area if the secondary structure 
stays intact. There are additional questions on how much 
other factors, like, e.g., DNA-pore interactions and elec- 
trostatics contribute to the blockade. 

To get an idea of the contribution of volume exclusion 
to the blockade current, we can use a st eady state form 
of the Nernst-Planck ( NP) equation (iBarcilon et al 



Schuss et all l200l[ ) 



19921: IChen e; Il992l: ICoalson and KurnikovaL 12005 : 



Jz = efinEz — eD 



dn 

dz 



(9) 



for the charge flux in the z-direction, where e is the 
electric charge (assuming monovalent cations), /i is the 
ion mobility, n is the charge carrier density, D is the 
diffusion coefficient, and Ez is the driving field in the 
z-direction. The general assumption that goes into the 
NP equation is that of a continuous ionic distribution 
contained within a continuum dielectric. The particular 
form of the Eq. [9] assumes that there is no free energy 
barrier to ion transport through the pore and thus it is 
a great start to examine whether volumetric effects are 
responsible for the blockade current. 

In the absence of any substantial diffusive term, the 
open pore current reduces to 



eiinEz 



(10) 



where Ez ~ V/Lp (the voltage V drops mostly over the 
pore length Lp). The quantity efin = cr is the conductiv- 
ity of the pore. 

If volume exclusion were the sole factor in the change 
in current, one could write. 



rib 



Vp-Vn 
— 1 



Fn 



(11) 



where rib stands for the carrier density during a blockade 
event, Vp is the pore volume for one repeat unit of the 
PN, and Vn is the volume of a nucleotide (see Table HI. 
For the fraction reduction in current from open (/q) to 
blocked (lb) pore we then get 



I- F . 



(12) 



More simply, if one assumes that the carrier density is the 
same, but that part of the pore area Ac has been blocked, 
then (/o — Ib)/Io = Ab/Ao, with Ab the portion of the 
area which is blocked. If we assume the same length I 
for both the unblocked and blocked pore, this relation 
will also give equation 1121 The whole analysis requires 
a regularity in the orientation of the blocking species, 
which may not be completely correct. 

For random coil (extended) ss-DNA, F can be matched 
with the volume fraction from Table [D The volume of 
the nucleotides in Table [J are for DNA nucleotides, but 
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should be essentially the same for RNA nucleotides, and 
the repeat units in the PN, whether charged or passi- 
vated. This is because i) the extra hydroxyl group has 
marginal effect on molecular volume, and ii) the sec- 
ondary structure should also not change the volume oc- 
cupied per nucleotide. The reason for the latter is simple: 
the molecular volumes were computed with the van der 
Waals radii, which are of the same order as half the base 
separation in stacked secondary structures. For example, 
the van der Waals radius is 1.7 A, but the base separation 
IS ^ 3 A. Thus, one may expect that secondary structure 
can slightly increase or decrease the nucleotide volume, 
but only by a small fraction. 

Consider now random coil poly(C), or poly(dC), within 
the a-hemolysin pore. Let us take 15 A as an esti- 
mate for the pore radius, and suppose the PN is max- 
imally linearly extended such that each nucleotide occu- 
pies a 7 A length of pore.^ In this case, I — F = 0.26 
for cytosine, whi ch does not come cl ose to the blockade 
values (> 90% (jAkeson et all Il999l )'). When the sec- 
ondary structure changes, this will change the length of 
the pore that each nucleotide occupies. If, for instance, 
the bases stack a t a dis tance 3.4 A (poly(C) stacks at 3.1 
A (|Arnott et all I1976D ). giving a minimum occupation 
length in the pore, then 1 — F = 0.52. This is still much 
less than the blockade values found experimentally. 

Of course, one may need to compute molec- 
ular volu mes including the hydrat ion of the nu- 
cleotides ( Deamer and Brantonl |2002| ) and also divide 
the pore into "good" and "bad" volume, as there may 
be bound water molecules on the walls. This would cer- 
tainly give a larger fraction reduction in volume. How- 
ever, defining volumes with the hydration layers included 
may not be accurate in the dynamical pore environment, 
as translocating ions may temporarily share hydration 
layers with the nucleotides. But this raises some interest- 
ing questions about the interaction of the different bases 
with the nearby water and ions. For instance, do the 
bases interact sufficiently different with the water that 
the free volume, including hydration layers, distinguishes 
the bases better than in Table HI or is it the reverse? 
One would expect the latter, because the hydration lay- 
ers probably smooth over differences in the bases, and 
the backbone is likely the main location for bound water 
molecules. 

Regardless, the argument seems strong enough that it 
points to other effects (interactions with the pore walls 
and electrostatics) to be the main cause of the block- 
ade current for DNA. Other polymers, especially neu- 
tral ones such as polyethylene glycol, may show differ- 
ent behavior. However, the the blockade is still a com- 
plex phenomenon. For example, poly(C) and poly(dC) 
have the same base, but the interaction with the pore 



walls and dynamics could change due to the secondary 
structure. In this case, the more structured poly(C) 
may create dynamical charge traps within the pore (a 
steric effect), thereby reducing the current further than 
for poly(dC). Also, poly(dC) creates larger current block- 
ages than poly(A) (89% and 86% for poly(dC) co mpared 
to 85% and 55% for poly (A) (jAkeson et a/.l . ll999t) ). Since 
C has a much la rger asso c iation with the lysine in the a- 
hemolysin pore ( Bruskovl . ll978l ). interactions may be the 
determining factor. 

Further, there are other effects to consider. Volume 
exclusion creates an additional electrostatic barrier due 
to the increased confinement of the ions. Thus a small 
change in volume can decrease the current because the 
associated free energy change "deactivates" many ions for 
transport (see bel ow). Addit i onally, PNs remain charged 
within the pore (|Cuil . 120071 : iKevser et all . l2006l) . How 
does this affect the blockage of colons, and the flow of 
counterions?"'" We are not aware of a complete answer to 
this question. 



2. Electrostatics 

Another issue related to the ionic currents is how the 
electrostatic environment of the pore affects transport? 
Some partial answers can be found by looking at simpli- 
fied models of the pore environment, and using this as 
input to the full NP equation.^ For instance, one could 
take an appropriate free energy potential, and possibly 
position-dependent dielectric constants and mobilities. 
In the case where the pore diameter is small, the free 
energy change can be quite large and thus significantly 
suppress ion transport. 

In the regime of large free energy change, we consider 
a scenario where the two chambers have an equilibrium 
concentration of ions at their respective shifted chemical 
potentials. Most of the ions do not have the energy to 
transport across the pore. So we consider an activated 
process where the number of ions available for transport 
is given by 



n — 



'AFb/kT 



(13) 



where Uq is the bulk density and AFi, represents the free 
energy barrier due to both an electrostatic energy barrier 
and an entropic barrier. For large pores, both of these 
contributions will be negligible as the electric field from 



t This would give about 7 nucleotides in the pore stem, but there 
may be more depending on secondary structure. 



t An unusual phenomenon of current enhanc ement due to DNA 
transl ocation wa s found in some experiments llChang et al. 1. 12004 
I2OO6I: [Fan et all 2005). This has been explained in terms of an 
enhanced counterion current due to the presence of the backbone 
charge of the DNA, which at low ionic concentrations overpowers 
the volume exclusion. 

§ One can also question the validity of that equation because in- 
homogeneities at the nanoscale may be important. We are not 
aware of a critical study of this point. 
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FIG. 15 Continuum pore model used by IZhang et all |2005| . 
The dielectric constants are such that ki ^ «2 so that all field 
lines of a charge in the pore run down the pore axis. The lower 
panel shows the electrostatic energy barrier for a single charge 
to transport from the bulk into the pore. From lZhang et aU . 
I2b05i . 



the ion is still embedded in a large volume of water and 
there is still a large phase space for ions to occupy. 

For smaller pores, one can think of the simplified 
model shown in Figure 1151 The general concept is 
to look at the pore as a continuum large dielectric 
(water) region of space surrounded by a material of 
low dielectric (Si02,Si3N4), thus creatin g an e l ectro- 
statically qu a .si-one -dimensional structure (|Tebeij . l2005t 
IZhang et all l2005f ). Under this simplified picture, the 
electrostatic and entropi c components to AFi, can be 
computed. IZhang et all [2005 found that this simpli- 
fied model gives a tremendously suppressed ionic cur- 
rent due to the large free energy change of ions enter- 
ing the pore. However, they also found that any sig- 
nificant presence of wall charges, which are present in 
experimental systems such as silicon dioxide, would de- 
crease the fre e energy barr i er and thus increase the ionic 
conductance ( Zhang et all |2005|) . Many recent papers 



this 


point of view ( 


Bonthuis et all 20061: Kamenev et al. 


2006 


; IZhang et al. 


. 20051 2006al: IZhang and Shklovskii 



20071) . 

In our own work ( Lagergvist et all l2007c|) . we have 
looked at the microscopic effects of the electric fields and 
solution structure. If one examines the field of a single 
anion surrounded by a spherical droplet of water, this 
looks very different than one would expect from a bulk 
dielectric.^ The field obtained from MD simulations, and 
averaged over time and solid angle, is shown in FigurefTBl 
One can clearly see the formation of a nanoscale struc- 
ture around the ion: the water molecules start to form 
layers (called hydration layers). For instance, for the Cl^ 



^ This is, of course, due in part to the fact that macroscopic elec- 
trostatics requires averagin g the electric fi eld over large regions 
of space of at least 10 nm llJacksonl . [19981) . It is also due to the 
formation of the hydration layers as a consequence of the strong 
local electric field. 




FIG. 16 (Color in online edition) Microscopic electric field 
of a single anion in a spherical water droplet of radius 25 A 
(solid black line). The dashed blue line is the bare field. The 
inset shows a model of the first hydration layer. 



we show here, there are on average six water molecules 
in the first layer which orient their dipoles very strongly 
toward the anion. These hydr ation layers only get par - 
tially destroyed within a pore ( Lagergvist et a/.l . l2007cl) . 
and the free energy barrier has an electrostatic contri- 
bution given by the partial destruction of the hydration 
layers. As the pore radius is increased, there is a critical 
radius that allows a full hydration layer to be transported 
through the channel. This creates a strong non-linear de- 
pendence of the free energy barrier as a function of pore 
radius and also ionic concentration. Effects like the above 
illustrate the need to take a microscopic approach to un- 
derstanding nanopore electronics, and foreshadows novel 
physical phenomena that will be observed at the inter- 
face between solids and liquids (and biomolecules) in the 
nanometer regime. 



V. SEQUENCING AND DETECTION 

We have reviewed the physical characteristics of the 
bases and polynucleotides, and also discussed nanopores 
and DNA translocation. Thus, we are now in a position 
to examine in more detail the physical mechanisms by 
which DNA can be detected and sequenced. Very loosely, 
we divide the physical mechanisms into three categories: 
electronic, optics, and force methods. In this section we 
will discuss several of the proposed methods in terms 
of the physical differences of the bases, and outline any 
experimental and theoretical results pertinent to either 
measuring these property differences or constructing a 
working apparatus. 
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(see H asianowicz et al\. Il996l iDeamer and Brantonl . 
I2OO2I and iDeamer and Akesonl ! |2000D . With this idea 



FIG. 17 (Color in online edition) (a) Engineered a-hemolysin 
pore used to detect nucleotides via their ionic blockade. 
Adapted from lAstier et al. I I2OO6I . (b) Schematic of two elec- 
trodes embedded in a nanopore. As PN translocates through 
the pore, the electrodes drive a current across the nucleotides 
in the transverse direction. The current across each of the dif- 
ferent nucleotides provides an electronic signature of the base. 
E|| is due to the voltage pulling the DNA through the pore. 
Ex is the field perpendicular to the electrode surface and is 
provided by the applied voltage across the electrodes and any 
additional external capacitor. From Lagcrqvist et al., 20073. 
(c) Middle panel: schematic of a nanopore through a capacitor 
made of doped polycrystalline silicon-Si02- doped crystalline 
silicon . The Si02 can be as thin as 0.7 nm (|Gracheva et all 
l2006bl ). which is ab out the spacing of nucleotides in extended 
ss-DNA. From .Gracheva et aUr2006bl . 



A. Electronic 



The electronic methods proposed so far 

are based on the ionic blockade c urrent in 

the nanopore ([Deamer and Brantonl . 120021 : 

iKasianowicz et al. I, Il99d) . embedding nanoscale elec- 
trodes in the n anopore to measu r e transverse transport 
acros s ss-DNA (iLagerqvist et al. I J2006I : E ee and Melleil . 
120071: IZwolak and Di Ventral . 12005'), or measuring the 
voltage fluctuations in a cap acitor across th e longi- 
tudinal direction of the pore ( Gracheva et all . l2006bl : 



iHeng et all l2005al) . Three specific proposals using these 
techniques are shown in Fig. [T71 



1. Ionic blockade 



Wi th the promising first results of IKasianowicz et all 
Il996l . it was speculated that it may be possible to 
sequence DNA by measuring the ionic blockade in 
a very particular pore under the right conditions 



in mind, researchers started to examine what kind of 
information can be extracted from the nanopore exper- 
iments. Indeed, other experiments with a— hemolysin 
pores have shown that the ionic blockade can be used 
to detect blocks of nucleotides, and that the blockade 
characteristics give information on secondary structure 
and directionality (5', 3') of PN translocation. 

For instance, by looking at RNA hom opolymers 
poly(A), poly(U), and poly(C), lAkeson et a/.l .il999 have 
demonstrated that the blockade current and the translo- 
cation duration can be used to distinguish between the 
type of base present in the homopolymers. At 120 mV 
bias, it was found that poly(U) translocation takes 1.4 
fis/nt and 6 /its/nt for two translocation event types com- 
pared to 22 /xs/nt for poly(A). Yet, they both can give 
blockade currents of essentially identical magnitude (ap- 
proximately 85% blocked). On the other hand, poly(C) 
was found to give 95% and 91% blockade of the ionic 
current, and to translocate at ~ 5 /is/nt. 

An interesting feature of the above results is that 
poly(C) gives a larger blockade of the ionic current than 
poly(A) even though it is a smaller base, see Table |T| 
lAkeson ~et "aLl. ir999 have suggested that this could be due 
to the secondary structure of poly(C), which at neutral 
pH and room temper ature ha s a hel ical structure that is 
1.3 nm in diameter (|Saengej . Il988l) . This size is small 
enough to fit into the a-hemolysin pore without unrav- 
eling. To test this, the above authors have compared 
poly(C) to poly(dC), its DNA counterpart. Poly(dC) 
secondary structure is not as stable as poly(C), so differ- 
ing results for poly(dC) and poly(C) would give evidence 
for secondary structure as the cause. lAkeson et al. I. I1999I 
indeed found that poly(C) blocked more ionic current 
than poly(dC). Thus the secondary structure of poly(C) 
is the likely cause of the larger blockade. This, however, 
leaves an open question raised earlier in section IIV.DI 
why does the secondary structure cause a larger block- 
age of current? 

The slow translocation of poly(A) is also thought to 
be due to secondary structure. Poly(A) has a sec- 
ondary hehcal structure, that, unh ke poly(C), is too bi g 
to translocate through the pore (|Akeson et all Il999l ). 
Thus, poly(A) takes extra time to unravel and go 
through the pore. On the contrary, poly(U) does not 
have any secondary structure under the reported exper- 
ime ntal conditions (see discu ssion in lAkeson et al. I. I1999I 
and IKasianowicz et aLl . ll996D . 

It has also been shown that a PN with two homoge- 
neous blocks, A30C70, gives ionic blockade events with 
a stair structure: a higher current is present when 
the homogeneous block of A's is present in the pore, 
but t his current d e crease s when the C's are in the 
pore ( Akeson et all Il999l ). Although far from single- 
nucleotide resolution, this shows that one can obtain 
some internal information about the strand from the 
nanopore. This result together with using the aver- 
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age blockade current, av erage translocation time, and 
the temporal dispersion ( Meller et al\ . |2000| ) to distin- 
guish individual molecules, shows some of the promise of 
nanopore technology in polynucleotide detection. 

Later, gea et al ]. l20n5al took advantage of the abil- 
ity of solid-state pores to work under many conditions by 
examining DNA translocation at many different pH's in 
order to detect the ds-DNA denaturation into ss-DNA. 
As they increase the pH from 7 to 13, these authors find 
that there is a rather abrupt transition when the current 
blockade drops by roughly a factor of 2 (Fologea et all . 
l2005aD . That this is due to denaturation is confirmed 
by measurement of the optical absorbance which shows 
that denaturation occurs at around pH=11.6. These 
are examples of physical processes that can be studied 
with, and detection capabilities that can be achieved by, 
nanopore technology. 

In order to actually sequence a strand of DNA single- 
base resolution is required. In this case, the different sizes 
of the bases, and the different interactions between the 
bases and pore, have to be detected. For the direct use 
of the ionic current as an indicator of the base present at 
a location in a polynucleotide (i.e., for sequencing), the 
nanopore must have a length Lp of roughly one nucleotide 
(< 1 nm) and diameter somewhere in between 1 nm and 
2 nm. 

If the ionic blockade is solely due to excluded volume, 
then one would expect that the differences of the bases in 
a 1.5 nm diameter, 0.7 nm long pore would be just a cou- 
ple of percent, as shown in Table HI The noise on the cur- 
rent itself, both intrinsically (due to ionic fluctuations) 
and due to structural fluctuations of the nucleotides is 
likely to be much larger than this.* For instance, the fluc- 
tuations of the ionic blockade current for homogeneous 
sequences is about 30% of the average current (see Figs, 
in . Mel ler et ai, 2000). The noise for a single base should 
be larger than for the sequences. Intuitively, comparing 
to an Qf-hemolysin pore with about 10 bases in the stem, 
one expects the noise for one base to be larger by a factor 
of about a/IO. This would mean noise on the order of a 
100% of the average current. There may be possibilities, 
such as changing conditions like temperature, pH, etc., 
to minimize secondary structure effects and make the cir- 
cumstances more amenable to sequencing. But even in 
the best case scenario, where only the pore-located nu- 
cleotide is controlling the ionic blockade current, it may 
not be possible to successfully distinguish the bases. See 
further the arguments presented in Sec. IIV.DI 

Nonetheless, more recent experiments have shown that 
single ba se differences can be de tected in strands. For 
instance, IVercoutere et all . 1200 ll have looked at several 



Molecular dynamics simulations of lAksimentiev et aZ.I . [20041 in- 
dicate that fluctuations in the ionic current due to structural 
changes are larger than the differences between bases. These 
structural changes are in part caused by interaction with the 
pore surface. 



different length DNA hairpins that contain no single- 
stranded portion. The hairpins initially block the 
vestibule part of the a-hemolysin pore, and would only 
be pulled through the pore when the double strand tem- 
porarily unravels. These authors have found that one can 
discern two hairpins with only a difference of one base in 
the loop and also two hairpins with only a mismatch as 
a difference. Both the vestibule blockade current and 
the vestibule blockade duration were found to be differ- 
ent. The hairpin with the base pair mismatch is much 
more likely to unravel, and therefore would spend much 
less time in the vestibule before translocation (a factor of 
100 less time according to these experiments). Similarly 
for the shorter loop hairpin, which is under much more 
strain. 

Using ionic blockade, the same group has looked 
further at hairpin differences, computer- learning algo- 
rithms for distinguishing signatures , and the physics 
of hairpins with i n the pore s IDeGuzman et all . 120061 : 
IVercoutere et all . l200ll l2003t IWinters-Hilt et all . l2003f) . 
Another group obtained similar findings for a different 
type of hairpin experiment, where a single base differ- 
ence in a single-strand leg o f a hairpin could be de- 
tected ( Ashkenasv et all . l2005t ) . It is important to note, 
however, that these experiments do not actually achieve 
single-base resolution as required for sequencing. 

Without amplification of the bases, it is unlikely that a 
strand of bases can be sequenced with the bare ionic cur- 
rent. However, there is potential to go beyond the nor- 
mal c urrent blockade experiments by creating des igner 
pores (|Bavlev and Cremeil . [200lHSiwv et aLl . l2005f ). For 
instance, introducing a foreign molecule into the pore or 
by using a pore that has some chemically-specific affin- 
ity, the pore can interact in a particular way to each of 
the four bases of DNA. A specific interaction could, for 
instance, create a different level of pore blockade for each 
of the bases. In this way, a distinguishable ionic current 
may be obtainable. The u nderlying premis e was demon- 
strated some years ago by iGu et all . Il999l with the use 
of an "adapter" molecule in a a— hemolysin pore. The 
adapter molecule was non-covalently inserted into a pore 
and it helped detect organic molecules that go through 
the pore by specific binding interactions which change 
the value of the current blockade and the translocation 
time. 

A very recent article has gone a long way in demon- 
stra ting the use of a, molecular adapter in sequenc- 
ing ( Astier et a/.l . l206^ . These authors have used a mu- 
tant a-hemolysin pore, see Fig. I17r a). with a positively 
charged cyclodextrin adapter. The adapter can bind and 
unbind from the pore (it is not covalently bound), cre- 
ating a stochastic signal of successive "on" and "off" 
current blockades. The addition of nucleotides changes 
this stochastic signal, creating additional events with a 
smaller blockade current when the nucleotide binds to 
the adapter. These smaller blockade currents are sig- 
nificantly different for the different bases, as shown in 
Fig. [THl Although the apparatus is chemical, it is really 
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FIG. 18 (Color in online edition) (A) Current trace with 
all nucleotides present in solution (time is in millisec- 
onds) ;_(B}_CuTrent_d^^ from the engineered pore. 
From lAstier et ail |2006| . 



detecting a physical difference between the nucleotides. 
The difference between the bases results in ionic current 
distributions with partial overlap, which gives the ability 
to correctly identify a base to about 93-98 % accuracy. 
Thus, multiple measurements of each base in a sequence 
are necessary. Unlike the transverse transport ideas 
we discuss below, multiple measurements may not be a 
controllable feature of the procedure and therefore the 
present technique may require re-sequencing each strand, 
or using some sort of oscillating voltage to pull the nu- 
cleotide in and out of the adapter to obtain the desired ac- 
curacy ('Astier et al'., '2006'; 'Pi Marzio and Kasianowic3 . 
[2003; Kasianowicz et al., 2002). 

A full sequencing procedure might employ exonucle- 
ase digestio n to take a DNA s trand and remove a base 
at a time ( Astier et all . \200(i) . The pore itself would 
have to have the exonuclease attached so that each base 
is released into the pore and its almost unique blockade 
signature detected. In addition, in order to actually cou- 
ple the exonuclease digestion with the pore detection, one 
has to deal with the fact that the adapter is not covalently 
attached to the pore and can unbind. Also unknown is 
whether the nucleotides actually pass through the pore or 
whether they just bind and unbind from the adapter stay- 
ing in the cis chamber (where they are added). Further, 
the exonuclease could create the rate-limiting step. For 
exampl e, a fast digestion rate of 1000 nt/s has been ob- 
served ( Matsuura et all . hoOlh . which is still much slower 
than the intrinsic speed of the nanopore experiment, on 
the order of 10^ nt/s. However, the rate will also depend 
on the binding timescales and the bandwidth of the ionic 
current measurement. 

It is clear from the experiments described in this 
and e arlier sections (a s well as results not dis- 
cussed ( Wang et a/.l . l200^ ). that the detection capabili- 



ties and the ability of ionic current experiments to probe 
physical processes at the nanoscale is a promising field. 
The use of molecular adapters to create a working se- 
quencing method is also promising, but there may need 
to be additional advancements to increase its intrinsic 
speed. 



2. Transverse electronic current 

When measuring the ionic current the size, secondary 
structure, and base-pore interactions distinguish the 
bases. One can also measure intrinsic electronic prop- 
erties of the bases by, for instance, embedding elec- 
trodes within a nanopore to measure the transverse 
curre nt through ss-DNA as it tran slocates through the 
pore (jZwolak and Pi Ventral l2005l ). This method es- 
sentially detects the electronic structure of the bases 
(i.e., the degree of delocalization, shape, and energies of 
the bases' electronic states) relative to the specific elec- 
trodes. A schematic of the idea is shown in Fig. [TTTb). 
We stress that transverse transport is different than the 
(also interesting and yet not f ully solved) problem of 
longi t udinal transport in DNA (iDi Ventra and Zwolakl 
l2004l : lEndres et a/I 12004 iPorath et all |2004) . The 
latter is another complex problem involving a disor- 
dered quasi-one-dimensional molecule. Further, mea- 
suring transverse currents in other scenarios has been 
suggeste d as a way to explore base differences and ne w 
physics (|Apalkov and Chakrabortvi l2005t iMacial . |2005[) . 

Initial calculations on individual nucleotides ide- 
ally configured between two nanoscale gold electrodes^ 
(see Fig. IHb)), indicate that the different bases 
can g ive quite different cu r rents under some condi- 
tions (jZwolak and Pi Ventral . |2005| ). However, the dif- 
ferent currents are very sensitive to the electrode spac- 
ing, and also vary greatly if the nucleotides have com- 
plete freedom in orientation. For a 1.5 nm electrode 
spacing and the bases standing upright with respect 
to the electrodes (Fig. [5l^b)), the electronic currents 
for the nucleotides are different by orders of magni- 
tude, except for the nucleotides G and C which dif- 
fer only by about a factor of 2. The reason for this 
large d ifference is the coupling of th e bases to the elec- 
trodes (IZwolak and Di Ventra . I2005D . Due to the differ- 
ent sizes of the bases, as shown in Section HIIl their molec- 
ular orbitals have a different spatial extent. The largest 
base, adenine, thus couples better to the upper electrode. 
More recently the relation between base orientation/size 



t Note that the electrodes need not be made of gold. For in- 
stance, Golovchenko and co-workers are pursuing the same idea 
with nanotube electrodes. The nanotube electrodes are likely to 
be less reactive with the liquid and solid environment and may 
thus be more stable than metals commonly used for contacts. 
One may also use the n anotube as part of other detection appa- 
rati l lMeng et a^l2007^ . 
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and co nductance was shown explicitly in I Zhang et all . 
I2006bl . This is both an advantage, in that it gives added 
differentiation between the bases, and also a disadvan- 
tage, because any fluctuations in the orientation of the 
nucleotide can drastically change the current. 

To quantify this, let us look at the tunneling current 
through a single energy level En in the limit of zero bias 
V. For weak coupling, as is the case here, this is given 
approximately by 



TLiEF)TRiEF) 
nh {Ef - EmY 



V. 



(14) 



This equation shows that the current is proportional to 
the two coupling strengths Fl, Fij, and inversely pro- 
portional to the distance between the energy of the state 
and the Fermi level Ef- The different bases have nearly 
the same energy levels compared to the Fermi level (see 
Fig. [5]), thus this portion is not likely to provide much 
distinguishability. However, the coupling elements can 
be very different for the bases because of their different 
spatial extensions and wavefunctions (see Sec. IIILB]) . 

A preliminary look at the effects of changing orienta- 
tion showed that the transverse currents are still quite 
different exc ept for G and C, thus bolster ing the case for 
sequencing (jZwolak and Di Ventr i 120051 ). Also, unlike 
the voltage fluctuation measurements discussed below, 
the transverse current is minimally affected by nearest- 
neighbor bases because the current is controlled by base- 
electrode coupling and the energy of the molecular states. 
These are only slightly influenced by the neighboring 
bases, so long as the electrode width is on the order of 
the base spacing (~ tA for extended ss-DNA) . The base- 
electrode coupling is particularly important, and this 
drops exponentially with the distance of the bases from 
the electrodes. This is unlike the capacitance method 
we discuss below where the Coulomb interaction is long 
range''' and thus nearest neighbors can influence the volt- 
age signal. However, at such an electrode spacing of 1.5 
nm the current is already very small, on the order of 
picoamps for the base A. It is thus necessary to go to 
smaller electrode spacing and also examine the full effect 
of structural fluctuations. 

To obtain a more realistic look at the effect of 
structural fluctuations the two authors together with 
J. Lagerqvist have performed molecular dynamics 
(MP) simulations of ss-DNA being pulled through the 
pore (jLagergvist et aLl . l2006t ). In the absence of any con- 
trol on the DNA, the orientation of the bases varies wildly 
(the bases can be at any orientation, perpendicular or 
parallel to the electrode surfaces). This causes orders of 
magnitude fluctuations in the value of the current. Due 



-f At these atomic scales, and in particular inside the pore, screen- 
ing by the few water molecules present around the bases is not 
as effective as one would o btain by assumirig the m acroscopic 
dielectric constant of water jLagergvist et alWmf?i ) 



to these large fluctuations and small values of the cur- 
rent magnitude, it is thus unlikely that sequencing will 
be possible in the absence of any control on the DNA 
translocation. 

However, MD simulations show that a relatively weak 
transverse field can orient the nucleotide in the junction 
with respect to the electrodes in less time (in the lOO's of 
picoseconds) than it ta kes the nucleotide to translocat e 
through the junction (jLagergvist et l200fil liooii). 
Referring to Fig. [TTT b). the condition on the transverse 
field to be strong enough is E|| <C Ej^. This ensures a 
slow enough DNA velocity to allow the transverse field 
to orient the nucleotides while they are facing the elec- 
trodes. The transverse field strength can be of the same 
order of magnitude as that driving the electronic current, 
but may be provided by an external capacitor around the 
nanopore device. 

When the control is exerted on the nucleotides, one 
obtains current distributions for the different bases that 
are sufficiently disjoint, see Fig. [TH so as to allow the 
bases to be statistically distinguishable. In this case, the 
statistics can be gathered both intrinsically and extrin- 
sically. First, the strand has some finite velocity, thus 
allowing for multiple measurements to be gathered as 
the nucleotide translocates through the junction region 
(so long as the inverse of the bandwidth of the electrode 
probes is much smaller than the translocation time of 
a nucleotide). Second, and more importantly, the finite 
bandwidth of the probes itself samples over many config- 
urations of the intervening nucleotide. 

Overall this leads to a sequencing protocol where one 
first has to measure the current distributions for each of 
the nucleotides using a homogeneous strand. These dis- 
tributions will depend on the particular pore geometry 
and conditions and thus will be particular to each device. 
Then one can sequence the DNA by pulling it at a ve- 
locity that allows each nucleotide to stay in the electrode 
region for enough time to collect a current distribution. 
Then by comparing with th e "target" distributions , one 
can determine the sequence ( Lagerqvist et a/.l . l2006l ). To 
get an idea about the potential speeds achievable with 
this method, we estimated that a single run through the 
3 billion bases on a single-strand of human DNA would 
take a raw time of about 7 hours, without paralleliza- 
tion ( Lagerqvist et aLl . l2006[ ). 

There are other issues yet to be explored, some of 
which have to be addressed experimentally. For instance, 
the construction of a nanopore with embedded electrodes 
remains a formidable challenge to the implementation 
and testing of the method described above. However, 
proof of principle of the distinguishability of the bases via 
transverse current measurements can be done with other 
methods, such as using a scanning tunneling m icroscope . 
This r oute has been explored in the work of IXu et alV 
l2007allbl . These authors have shown experimentally that 
the four bases (using isolated nucleotides) provide a dis- 
tinguishable signal via their HOMO level when mea- 
sured with an STM. Much earlier, other groups have also 
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FIG. 19 Current distributions for the different nucleotides in 
a 1.25 nm electrode spaced pore and at 1 V, with a transverse 
controlling field. The inset shows the error of misidentifying 
a base versus the number of independ ent counts (or measure- 
ments). From iLagergvist" et all |2006| . 



examined isolated bases or bases in a strand of DNA 



with scanning probe techniques (Hamai et all 200C ; 


Tanaka et al.. 19991: Tanaka and Kawai 


2003; 


Tao et al.. 


19931 Wang et aL.boOl). The idea of 


using transverse 



transport to sequence is being p ursued ex perimentally by 
a few groups (M. Ramsey et al. ('NHGRl', |2006l). Bran ton 
and Golov chenko ([Branton arid Golovchenko . fWEBl ) . J. 
Lee et al. (|Lee and Melleil . |2007| )T 

There are, of course, drawbacks to the transverse trans- 
port approach we have described. Clearly, if the amount 
of noise is large, it will make the signals from the bases 
potentially indistinguishable. There will be 1/f noise 
and noise provided by ionic fluctuations (discussed be- 
low). While one can reduce 1/f noise by operating away 
from zero frequency, the role of ionic noise is less clear. 
Also, one may wonder in a real implementation whether 
clogging of the pore can occur, i.e., whether the small size 
of the pore and the transverse electric field(s) will cause 
too much interaction of the DNA with the surface of the 
pore. This is an issue that has yet to be investigated: its 
study requires a fully quantum mechanical treatment of 
the DNA motion in the pore. 

One very important and interesting issue, both funda- 
mentally and for the realization of the above approach, 
is the effect of ionic fluctuations on electronic trans- 
port. It has importance for many cases of charge trans- 
port i n soft materials, such as long i tudinal transport in 



port i n soit materials, such as long i tudinal transport in 
DNA (iDi Ventra and Zwolakl 12004 lEndres et all 12004 



iPorath et a/.l . l2004 ). It would also be interesting to inves- 
tigate what happens to the tunneling current in the pres- 
ence of fluctuating energy landscapes created by classical 
objects (e.g., ions) that cross the electrons' path, when 
the rate of these fluctuations is comparable to the inverse 
of the coherence time. An interesting experiment in this 
direction may be to use different weight counterions to 



observe the transition from slower to faster fluctuations. 
The effect of ions on transport is part of the more general 
problem of understanding the ef fect of changing a sin- 
gle atom in mole cular junctions I Di Ventra et all l2000l : 
|^n£liLaI],[200i). 

Finally, another geometry has been proposed to probe 
electrical transport perpendicular to the base planes 
(i.e., w ith the base plane pa rallel to the electrode sur- 
faces) (|Lee and Melleil . l2007t ). However, there are two 
challenging issues with this setup. For one, since the 
transport is now envisioned across the base planes, the 
electrode spacing has to be very small, < SA, to obtain 
a measurable tunneling current through the bases. On 
the other hand, this will make it difficult for the ss-DNA 
to translocate through the pore. Fro m our MP simula- 
tions with s elected initial conditions ( Lagergvist et all . 
l2006l l2007aD . the smallest pore that will allow DNA 
to translocate is lOA in diameter. In addition, even 
with an electrode spacing that allows for a measur- 
able current, the differences between the bases in the 
above planar conflguration are not as large in mag- 
nitude as in the case in which the bases stand up- 
right with re spect to the electrode surfaces (as shown 
in Fig. Eljb)) (|Zwolak and Di Ventral . [200l |2005^ . This 
is due to the fact that the bases have very similar HOMO 
and LUMO charge distributions when viewed perpendic- 
ular to their plane, whereas when standing upright, they 
have different couplings (see Fig. . 



3. Capacitance 

The measurement of voltage fluctuations in a metal- 
oxide-silicon capacitor combined with a nanopore has 
been proposed as a method to detect and obtain 
the length of DNA, and potentially to se quence 
it (|Gracheva et all l2006bl : iHeng et all l2005al ). A 
schematic of a nanopore capacitor is shown in Fig. ll7r c). 
As each nucleotide passes through the pore, the charge 
on its backbone can induce a voltage across a capacitor 
in the longitudinal direction. By measuring these volt- 
age fluctuations, one can effectively count the number of 
nucleotides. For such a scheme to be used for sequenc- 
ing it would have to also be sensitive enough to detect 
the dipole moment differences of the nucleotides given in 
S ection [ml 

iHeng et all l2005al have indeed observed voltage sig- 
nals on the two doped-silicon electrodes as a strand of 
DNA passes through the pore, the difference of which 
gives oscillations. The oscillations have a magnitude that 
could be adequate to detect the dipole difference of Ade- 
nine and Thymine (jHeng et al. I. l2005all . However, exper- 
imentally these authors were not able to resolve even the 
charge on the individual nucleotides because of the band- 
width and RC time constant of the probe, and the large 
size of the pore. 

The same group also performed simula- 
tions (Gracheva et all l2006alrbr). to calculate the 
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expected voltage fluctuations on the capacitor as a 
ss-DNA is pulled through a 1 nm radius pore. In these 
simulations, the DNA strand was not allowed to have 
conformational fluctuations (s imilar in style to th e 
calculation shown in Fig. 1 of iLagergvist et all . |2006| ). 
An example signal obtained from this simulation is 
shown in Fig. 1201 Interestingly, the maximum signal 
obtained is ~ 35 mV for the nucleotide, 30 mV for the 
backbone, and 8 mV for the base. In addition, as a 
nucleotide moves through the pore the corresponding 
voltage signal has been found to be influenced by up 
to three nearby nucleotides. Even though the bases 
themselves show different voltage signals, the larger 
backbone signal seems to dominate and give almost 
id entical signals for the d ifferent nucleotides (see Fig. 3 
in lGracheva et a^ll2006a[) . 

Another important factor in the capacitance approach 
is the dimension of the pore. Estimates suggest it should 
have a bout 1 nm radius in order to obtain an adequate 
signal ( Gracheva et aLl . l200 6a,bV The 1 nm radius of the 
pore would both force the DNA to be stretched as it goes 
through the pore, thus aligning the bases perpendicular 
the electrodes (which produces a better signal than if the 
plane of the bases were parallel to the electrode surfaces) 
and maximizing their distan ce from each other (for in- 
stance, for the sequences in iGracheva et al. I, l2006al the 
bases are about 7 A apart, double their distance than in 
ds-DNA or in ss-DNA with strong secondary structure). 
Also, the small pore excludes a lot of water, thus reducing 
screening. 

Theoretically, though, one wonders how much, in prin- 
ciple, the bases can be distinguished via their dipole mo- 
ments. From Section lllli the dipole moments of C and 
G isolated bases differ by just a percent (the differences 
between A and T, and those with G,C, are much larger). 
This, however, ignores the effect of the backbone dipole 
and the ion-counterion dipole, which are both very large. 
In addition, the large fluctuations of ions at the entrance 
of the pore may contribute substantially to the voltage 
fluctuations. Indeed, most of the voltage drop occurs at 
the location of the (high-resistance) pore. Thus, ions of 
opposite charge would build up on the two open sides of 
the pore entrance. When DNA enters the pore, it creates 
an even higher-resistance pore. As a consequence, there 
should be a fluctuation in the ionic concentrations near 
the pore ends, thus creating a temporary voltage fluctu- 
ation on the capacitor. The consequence of this effect on 
sequencing, together with conformational fluctuations of 
the DNA bases inside the pore, has not been studied yet. 

To get some perspective on how meaningful the dipole 
moments are, let us consider the dipole moments of other 
species present in an actual experi ment. Water h a s a ga s 
phase dipole moment of ~ 1.8 D (jClough et al\ . Il973l ). 
However, in the conde nsed phase, the dipole moment o f 
water is larger, ~ 3 D (|Silvestrelli and ParrinelloLfl999h . 
This is easy to understand: the hydrogen bonding that 
takes place between water molecules induces more polar- 
ization in the electronic distribution, causing an increase 
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FIG. 20 (Color in online edition) An example voltage sig- 
nal obtained from MD simulations/electrostatic calculations 
of a 11 base single-strand of DNA translocating through a 
nanopore capacitor. The dips correspond to the nucleotides. 
Experimental signals have been obtained of the voltage fluctu- 
ation s, but without the fine resolution of these simulated sig- 
nals ([Gracheva et all . l2006bl ) . From .Gracheva et al., ,2006a . 



in the dipole moment. This is actually true for other 
species as well: the interactions between the polar sol- 
vent and the solute (the FN) will induce larger dipole 
moments in the solute. More speciflcally, differences in 
how water interacts with the different bases could signif- 
icantly change their dipole moment (this is an open and 
important question). 

In addition, the ions present in solution contribute to 
the dipole moment fluctuations. Consider, for instance, 
the counterion on the DNA backbone. The values of 
the deoxyribonucleic acid dipoles given in Table |TT] are 
for a passivated backbone. The passivation is simply 
the addition of a hydrogen atom on one of the partially 
charged oxygen atoms on the phosphate group. This ad- 
ditional hydrogen has a bond length of about lA. There- 
fore, if we consider the effect of a bound counterion, 
we have to add an additional dipole moment of about 
1 eA , because the ion-oxygen bond would be about 1 
A longer. This additional dipole thus contributes 4.8 D. 
Further, molecular dynamics simulations indicate that 
the counterions on the DNA backbone fluctuate quite a 
lot ( Lagergvist et aZ.l . l2007al ). Thus, there could be noise 
larger in magnitude than the nucleotide or base dipole 
moments themselves, and much larger than the differ- 
ences between them. However, how this noise interferes 
with the signal from the bases is not clear, since at these 
length scales the fluctuations may or may not give a small 
average noise. 
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FIG. 21 (Color in online edition) A double strand ol DNA 
with magnified sequence is pulled through a nanopore. As 
the strand goes through, it unzips, refeasing the ffuorescently 
tagged single strand, which is then read optically. Courtesy 
of A. Meller. 



B. Optics 

Another proposed method for sequencing relies on op- 
tical methods and ds-DNA unzipping via the interaction 
with a nanopore. The idea is to tear a part "magnified" 
DNA and then read an optical signal ( Lee and Melleil . 
120071 ). This process is shown in Fig. 1211 

In more detail, one uses a biochemical process to create 
"designer" ds-DNA where each base is represented by a 
unique iV-base sequence (with N k, 20). These iV-base 
sequences are just two blocks of 7V/2 base homogeneous 
strands. The 7V/2 base strands are each labeled with a 
fluorescent tag on one end and a quenching molecule at 
the other, such that within the double strand each tag is 
always paired to a quencher. This prevents fluorescence 
when the DNA is away from the pore. When the ds-DNA 
is pulled into, for instance, an a-hemolysin pore, the dou- 
ble strand has to unzip into single st rands ( Mathe et al\ . 
I2006L 12004 [Sauer-Budge tt aZ.l . l2003D . This pulls the flu- 
orescent tag away from its quencher and allows the tag, 
and thus the base in the original strand, to be detected 
via optical means. 

Despite the necessity to work with "amplified" DNA, 
this method has tremendous potential for sequencing, 
and can be parallelized quite easily by creatin g several 
pores in the same device (jLee and Melled . 120071) . P relim- 
inary estimates give a raw sequencing rate of '--^ 1 — 10 
million bases/s for a parallelized nanopore detector. 



C. Force 

One important technique that may be useful in 
DNA sequencing and detection, and also single-molecule 
and nanopore studies, i s the use of optical tweez- 





Biotin 




IFf LdM 



FIG. 22 (Color in online edition) Optical t weezer used in 
conjunction with a nanopore. From lKevser et a/.„,2006. . 



mented setup using optical twee zers in conjunction with 
a nanopore is shown in Fig. [22] (iKevser et all |2006[ ) . 

In this experiments, A-DNA is attached to polystyrene 
beads via a streptavidin-biotin interaction (this is non- 
covalent pro tein-ligand association and in this case is 
very strong (|Di Ventra and Zwolakl l2004t) ). The dielec- 
tric bead can be tra pped at the focal po i nt of a tightly 
focused laser beam ( Neuman and Blockl 1200^^ . If the 
bead is pulled away from the focal point, a restoring 
force is exerted on the bead. For small displacements 
this force is linearly proportional to the distance of the 
bead from the focus, e.g., it resembles a spring obeying 
Hooke's law. Thus, by measuring the position of the 
bead from the laser focus, which can be done by mea- 
suring the reflected light off the bead, one can determine 
the net external force on the bead via 



(15) 



ers ( Neuman and Blockl . |2004[) . An already imple- 



where ktrap is the stiffness of the optical trap and 
AZ is the distance of th e bead from the focus of the 
beam (|Kevser et aLl . l2006l ). 

By bringing the bead close to a nanopore, one can then 
both control the motion of the DNA in the nanopore and 
measure the forces exerte d on it by the pulli ng voltage. 
By measuring these forces. iKevser et all . l2006l have deter- 
mined that the effective charge for each nucleotide pair 
on the double-stranded DNA is about 25% of its bare 
value, i.e., 0.5 e instead of 2 e. For 100 mV pulling volt- 
age, this effective charge gives a force on the DNA of 24 
pN. 

Clearly, the concept of coupling the nanopore setup 
with an optical trap would open up many possibilities 
to study physics and chemistry of single molecules. One 
can also imagine many possibilities for using this setup, 
or more complex setups (for instance with two optical 
traps, one on each side of the pore), for sequencing tech- 
nology and also to help in proof-of-principle prototypes. 
IKevser et ^ . 120061 have succeeded in slowing the translo- 
cation of ds-DNA down by five orders of magnitude com- 
pared to its free translocation rate - bringing it down 



23 



from ~ 8 mm/s to 30 nm/s. For ionic blockade sequenc- 
ing, the ability of the optical trap to slow and control 
the DNA translocation would allow for more definitive 
experimental measurements of the different ionic block- 
ade currents for the various bases, and may also allow for 
the stretching of DNA so that secondary structure effects 
may be reduced or removed. 

For transverse current and capacitance measurements, 
the controlled motion of DNA in the nanopore can also 
help in a number of ways. For one, it can slow down the 
DNA so that one can see the gaps between the bases and 
understand how to count the bases as they go by the pore. 
If one were to use a more complex setup, with two optical 
traps on each side of the pore, the pulling voltage could be 
completely turned off. This would enable the condition 
^ Ej_ necessary for successful transverse transport 
sequencing to be satisfied (see Sec. IV. A.^ . Furthermore, 
it would reduce ion flow through the pore to a minimum 
(only Brownian type motion would be present) and may 
thus reduce noise from ionic fluctuations for both the 
transverse current and capacitance. 

In addition to the measurement of the effective charge 
of DNA in solution, optical traps have provided the nec- 
essary control and precision force measurements for the 
study of DNA int eracting with RNA polyme rase at the 
single-base level ( Abbondanzieri et al . 2005^ and eve n 
single-molecule sequencing (|Greenleaf and Blockl[200i) . 



VI. CONCLUSIONS 

In this Colloquium we have examined a variety of dif- 
ferent proposals for accurate, rapid DNA/RNA sequenc- 
ing and detection. These methods take physical ap- 
proaches to detecting single bases in a sequence. Some 
of the basic ideas have already been implemented exper- 
imentally and shown to be useful for detection of strands 
of DNA and also as probes for "global" properties, such 
as length or homogeneous sequences. It is, however, still 
unclear as to whether these results can be extended to 
true single-molecule sequencing of DNA. 

We showed explicitly with ionic blockade and trans- 
verse transport that the physical differences of base prop- 
erties come in the form of differing signal distributions. 
This will hold true for all other approaches because of the 
electronic and structural similarity of the bases. Thus, 
a single, instantaneous measurement signal of a micro- 
scopically fluctuating nucleotide can not be expected to 
fall within a unique range. This suggests a statistical 
approach to physical DNA sequencing. In some cases, 
the statistics may be built into the measurement appara- 
tus via time averaging of the signals, as is the case with 
finite bandwidth electrical measurements. It is still an 
open question whether noise, in particular ionic noise, 
will drown out any signal difference between the bases, 
especially with capacitance or transverse transport mea- 
surements. Other noise has to be reduced in the inher- 
ently noisy nanopore (or other) environment, and the 



measurement device has to have sufficient resolution to 
see the differences in base signals. 

Another open issue involves achievable read-lengths 
and sequencing rates. Various lengths of D NA have been 
pulled through nanopores. For instance, IStorm et all 
l2005b| pulled through -100 kbp ds-DNA through a 10 nm 
pore. However, the authors do not know of a systematic 
study of possible read-lengths of ss-DNA being pulled 
through nanometer scale pores. Since the forces on DNA 
in the pore are small (on the order of lO's of piconewtons, 
much less than the forces necessary to break covalent 
bonds), practical matters like entanglement of ss-DNA 
outside the pore or nonlinear dependence of transloca- 
tion velocities on DNA length are likely to limit read- 
lengths. Additionally, the translocation rate of each DNA 
through the nanopore has to be smaller than the device 
bandwidth, allowing each base to stay long enough at its 
optimal position in the detection apparatus (e.g., in the 
active detection region, like between the two electrodes 
in a transverse transport approach). 

Progress in techniques and proof of principle is likely 
to proceed in p art via the use of au xiliary systems. Op- 
tical tweezers (iKevser et all 20061). a nano tube bound 
to an AFlvl jKing and Golovchenkol |2005[ ). and other 
ideas ( Krieger , |2006t ) could give rise to techniques to 
slow polynucleotide (or molecular) motion in the pore, 
and also offer the potential to get rid of the ionic solu- 
tion/electrophoretic pulling. This could be tremendously 
beneficial in demonstrating proof of concept and/or noise 
reduction. Further, we alluded to the distinction be- 
tween chemical and physical processes earlier. Chemi- 
cal processes change the atomic makeup or structure of 
a species. They potentially add slow extra steps to a 
sequencing procedure. However, they can add enhanced 
distinguishability as well. For instance, one could imag- 
ine using the amplification of the optical method to- 
gether with ionic blockade or transverse transport. Since 
the nanopores can detect differences in homogeneous se- 
quences, it may be possible to detect amplified or de- 
signer blocks of nucleotides. 

In addition to the technological motivation, there is 
also a scientific motivation for studying the detection 
step of a sequencing protocol. The process of directly 
detecting physical differences between the bases is fun- 
damentally different from the chemical/optical processes 
currently in use. Thus, not only does it represent a sci- 
entific challenge with the goal of low-cost and rapid se- 
quencing and detection, it is likely to increase our under- 
standing of biological and other molecules at the atomic 
scale, and may lead to off-shoot technologies. 

In this Colloquium, we have raised several important 
open questions. Now that fabrication techniques of solid 
state pores are maturing, experiments can be tailored 
to answer these questions more directly by systemati- 
cally changing pore dimensions, materials, surface prop- 
erties, and environmental (e.g., ionic) conditions. With 
adequately controlled experiments, nanopores could be 
used to probe inhomogeneities at the nanoscale, such as 
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fine structures in "dielectrics" (e.g., hydration layers) and 
other interesting physical phenomena. 

Because of the widespread technological and scientific 
interest, it is impossible to predict where this field will 
go from here. Undoubtedly, ingenious techniques and 
methodologies will be invented which extend our abili- 
ties to detect and sequence DNA well beyond what has 
been discussed. Some examples have been mentioned: 
integration of biological pores with nanoscale probes, use 
of molecular base discriminators with biological or syn- 
thetic pores, etc. The most promising fundamental re- 
search will sort through the many competing effects that 
exist in these complex systems, and increase our under- 
standing of physical processes at the interface between 
solids, liquids, and biomolecTiles down to the nanometer 
scale regime. We are confident that a physical approach 
to DNA detection will yield a wealth of information on 
physical and chemical processes at the nanoscale, and, 
hopefully, a rapid and low-cost sequencing technology. 
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